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SECTION  I 


INTRODUCTION 


Understanding  the  performance  of  high-power  lasers  requires 
a  model  of  the  system  response  to  external  perturbations.  In  the 
case  of  small  perturbations  about  some  equilibrium  value  the  sys¬ 
tem  response  can  be  characterized  by  a  (frequency-dependent)  trans¬ 
fer  function. 

As  the  simplest  example  of  such  a  transfer  function,  consider 
the  problem  of  relaxation  oscillations  in  a  conventional  laser.  The 
system  can  be  described  by  the  following  set  of  equations: 

if  =  c  (g  -  gc>*  (1) 


I?  =  -  S  U  +  —•)  +  P  (2) 

9t  t  4>s 


where 


=  system  flux 

g  =  system  gain/unit  length 

<p  =  saturation  flux 
s 

t  =  upper  state  deactivation  time 
c  =  speed  of  light 
gc  =  cavity  coupling 
P  =  external  pumping 

Let  us  now  assume  steady-state  pumping  with  a  small  per¬ 
turbation  ,  i .  e .  , 

P  =  P  +  P.  ,  (3) 

o  1 
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where  pQ  is  the  steady-state  component  and  Pj  is  the  perturbation. 
Similarly,  if 


4>  =  4>0  +  ^  i 


(4) 


g 


g  + 


(5) 


r  =  1  +  <$>  /4> 
o  Ys 


(6) 


we  have  the  following  equations  (after  Laplace  transforming,  de¬ 
noted  by  <Ms)  ,  etc.) 


In  this  case 


If  any  singularities  occur  with  Rea]  s  >  0,  the  system  has 
an  instability. 
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In  this  case,  the  roots  occur  at 


s 


=Vr-i: 


(10) 


i.e.,  at  a  frequency 


OJ 


o 


(r-1 ) 

T 


(ID 


with  a  damping  factor  of 


F  =  2r/i 


(12) 


and  no  instability  obtains. 

In  this  report,  we  will  show  that  a  realistic  transfer  func¬ 
tion  for  high-power  lasers  must  include  the  acoustic,  or  thermo¬ 
dynamic  degrees  of  freedom  as  well.  The  reason  for  this  inclusion 
can  be  seen  in  schematic  form  in  Figure  1. 

Gain  and  flux  perturbations  lead  to  density  perturbations 
via  local  heating.  This,  in  turn,  causes  refractive  index  fluc¬ 
tuations  and  flux  perturbations  which  close  the  cycle. 

We  will  show  below  three  separate  examples  of  this  expanded 
approach  to  the  system  transfer  function  and  system  stability.  The 
first  two  deal  with  C02  EDLs ,  although  the  formalism  can  be  gen¬ 
eralized  to  other  types  of  lasers.  The  third  deals  explicitly  with 
the  chemical  laser. 

We  shall  see  that  when  the  thermodynamic  degrees  of  freedom 
are  included,  the  possibility  of  a  system  instability  exists.  In 
fact,  such  an  instability  has  actually  been  found  experimentally 
in  the  case  of  a  cw  C02  electron-discharge  laser  (EDL) . 

We  shall  now  briefly  outline  the  remainder  of  this  report. 
Section  II  contains  a  discussion  of  an  instability  arising  from 
the  transfer  of  energy  from  an  initially  oscillating  cavity  mode 
to  a  different  transverse  mode  by  an  acoustic  wave.  The  approach 
to  this  problem  is  the  solution  of  a  perturbed  propagator  equa¬ 
tion,  i.e.,  an  integral  equation  for  the  change  in  the  electro¬ 
magnetic  wave  after  a  round-trip  traversal  of  the  acoustically 
perturbed  cavity.  Two  simplifying  assumptions  made  in  this 
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Figure  1.  Schematic  Indication  of  the  Coupling  of  Flux,  Gain, 
and  Thermodynamic  Degrees  of  Freedom 
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section  are  that  (a)  the  initially  oscillating  mode  has  a  uni¬ 
form  transverse  fluence  distr ibution  and  (b)  that  both  electro¬ 
magnetic  modes  have  the  same  natural  frequency  of  oscillation. 

It  is  shown  that  the  oscillatory  heating  which  occurs  can  result 
in  an  unstable  interaction  between  the  two  electromagnetic  waves 
and  the  acoustic  wave,  i.e.,  a  perturbation  whose  amplitude  in¬ 
creases  exponentially. 

In  Section  III  we  discuss  a  different  type  of  instability: 
alignment  sensitivity.  Variations  in  the  gas  density  along  the 
flow  direction  lead  to  both  gain  variations  and  slight  changes  in 
the  propagation  direction  of  a  light  mode.  We  show  that  this  can 
result  in  extreme  cavity  alignment  sensitivity. 

In  Semcion  IV  we  return  to  the  mode-medium  interaction 
problem  of  Section  II,  which  is  now  reformulated  in  the  language 
of  Brillouin  scattering  (differential  treatment).  We  show  the 
equivalence  of  the  differential  and  integral  formulations  of  the 
mode-medium  interaction.  We  focus,  however,  on  an  alternative 
instability  mechanism.  We  show  that  an  acoustic  wave  of  fre¬ 
quency  03^  tends  to  resonantly  couple  an  initially  oscillating 
mode  to  a  new  transverse  mode  whose  natural  frequency  is  off¬ 
set  by  from  the  natural  frequency  of  the  intial  mode.  This 
is  just  the  ordinary  Brillouin  effect.  We  derive  equations  of 
motion  for  the  amplitude  of  the  new  mode  and  for  the  acoustic 
wave.  The  secular  equation  for  the  system  is  derived.  It  is 
shown  to  be  unstable  in  the  limit  of  negligible  gain-flux  oscil¬ 
lations  (that  is,  with  a  ratio  of  cavity  flux  to  saturation  flux 
<<  1,  so  that  gain  saturation  may  be  ignored).  A  brief  discus¬ 
sion  is  given  of  the  nonlinear  saturation  of  this  process. 

An  idealization  of  the  mode-medium  interaction  theory  of 
Sections  II  and  IV  is  the  implicit  assumption  that  the  net  loss 
(=  loss  minus  gain)  of  two  cavity  modes  may  be  very  nearly  equal. 
(See  Appendix  A.)  At  the  beginning  of  this  program  it  was  realized 
that  no  self-consistent  calculation  of  the  eigenmodes  of  a  loaded 
(i.e.,  gain-saturated),  unstable  resonator  of  moderate  Fresnel 
number,  such  as  a  cw  electron-discharge  laser,  existed.  Such  a 
calculation,  which  clearly  is  quite  relevant  to  this  program,  is 
described  for  the  first  time  in  Section  V.  Using  a  new  computa¬ 
tional  technique  developed  at  AERL,  we  have  solved  for  the  self- 
consistent  eigenmodes  of  an  unstable  resonator  whose  gain  is 
saturated  by  the  presence  of  a  strong  oscillating  mode  (the 
"fundamental"  mode).  We  find  that  the  fluence  distributions  of 
the  various  eigenmodes  are  modified  by  the  saturation,  and  that 
the  eigenvalues  are  shifted  from  their  empty- resonator  values. 

In  particular,  the  empty-resonator  loss  degeneracy  which  exists 
at  certain  ha  1 f- in tegra 1  values  of  Fresnel  number  disappears  in 
the  presence  of  a  saturated  gain,  i.e.,  different  modes  are  found 
to  have  different  losses.  This  is  a  significant  new  result  which 
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suggests  that  it  is  not  advantageous  to  avoid  hal f- integral  values 
of  Fresnel  numbers  in  the  design  of  an  unstable  resonator  which 
will  be  operated  in  the  saturated-gain  regime.  Fi'om  the  point  of 
view  of  our  program,  the  results  of  this  section  suggest  that  the 
model  of  Sections  II  and  IV,  which  assumes  that  at  least  some  modes 
have  nearly  equal  net  loss,  is  somewhat  simplistic.  However,  the 
nonuniformity  of  the  gas  density  has  been  neglected  in  this  sec¬ 
tion;  it  is  possible  that  when  the  nonuniformity  is  taken  into 
account  the  net  losses  of  several  modes  may  be  very  nearly  equal. 

In  Section  VI  we  extend  our  instabi 1 i ty- theory  techniques 
to  a  relatively  simple  model  of  a  supersonic  chemical  laser 
(HF/DF).  The  model  corresponds  to  a  one-dimensicna 1  flowfield. 

The  fluid-mechanic,  chemical  and  optical  coupling  equations  are 
obtained  and  Laplace  transformed.  The  stability  of  the  system 
is  studied  by  use  of  the  Routh  Criterion,  which  allows  the  sign 
of  the  roots  of  a  polynomial  equation  to  be  determined.  Numerical 
calculations  for  typical  operating  conditions  of  the  model  laser 
show  that  it  is  stable.  The  important  result  of  this  model  is 
that  it  determines  the  transfer  function  of  the  system.  The 
transfer  function  can  be  used  to  predict  the  magnitude  of  the 
laser  response  to  any  external  perturbation. 

Finally,  in  Appendix  A,  wc  discuss  the  acoustically  per¬ 
turbed  optical  resonator  in  a  more  rigorous  fashion.  Using  the 
integral-operator  formalism,  wo  derive  tjic  equation  satisfied 
by  the  amplitudes  of  the  various  eigenmodcs  in  tho  presence  of 
the  acoustic  perturbation.  We  reformulate,  in  the  language  of 
quantum  mechanical  perturbation  theory,  the  difference  between 
the  instability  mechanisms  discussed  in  Sections  II  and  IV;  both 
follow  from  the  present  treatment.  We  suggest  an  alternative 
model  to  explain  the.  experimentally  observed  flux  modulation  in 
the  electron-discharge  laser:  one  in  which  only  a  single  eigen¬ 
mode  oscillates,  but  an  acoustic  standing  wave  both  modulates  the 
cavity  loss  and  is  driven  by  nonuniform  heating  caused  by  the 
nonuniformity  of  the  eigenmode  and  changes  of  the  cavity  flux. 
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SECTION  II 


A  THREE-WAVE  LASER  INSTABILITY 


A.  INTRODUCTION 

It  has  been  known  for  some  time^  that  high-power  CO  cw 
Electric  Discharge  Lasers  (EDL)  exhibit  output  instabilities. 

This  is  evidenced  via  a  spontaneous  oscillation  of  the  output 
flux  usually  resulting  in  full  modulation. 

An  example  of  such  behavior  for  a  typical  CO2  cw  EDL  is 
shown  in  Figure  2.  To  the  left  of  point  "a"  the  output  is  stable, 
while  to  the  right  of  point  "b"  the  instability  has  already  built 
up.  Eetween  "a"  and  "b"  the  output  may  be  considered  1  incar ly 
unstable ,  and  periodic  in  nature.  This  suggests  the  possible 
success  of  a  linearized  model,  with  the  (linear)  instability 
evidencing  itself  via  an  exponentially  growing  root  of  the  char¬ 
acteristic  equation  of  the  system. 

Corroboration  of  such  a  model  may  be  obtained  by  reviewing 
interf erograms  of  the  laser  taken  with  flux  extraction  (Figure  3). 
Also  shown  is  the  laser  output  flux  versus  time  with  the  points 
corresponding  to  the  three  interf erograms  as  indicated.  We  note 
the  presence  of  low  amplitude  density  (or  acoustic)  standing  waves 
between  anode  and  cathode  where  the  instability  in  the  output  flux 
is  already  evident.  It  is  as  if  the  laser  were  being  acoustically 
Q-switched. 

Further  evidence  of  an  acoustic  connection  can  bo  found  in 
the  period  of  the  output  oscillations.  Typical  oscillations  of  a 
gain-flux  system,  such  as  are  found  in  relaxation  of  a  gain- 
switched  spike,  are  of  the  order  of  a  few  microseconds.  As  can 
be  seen  from  Figures  2  and  3,  the  period  of  the  oscillations  is 
of  the  order  of  50  usee.  For  typical  anode-cathode  separation 
of  about  5  cm,  with  an  acoustic  mode  number  ~  5  (Cf .  Figure  3)  , 
the  acoustic  time  is  typically 


(5  cm/5) 

4 

3  x  10  cm/scc 


~  35  usee 
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Figure  2.  Evidence  of  Linear  Instability.  The  upper  trace  in  the 
top  figure  is  the  flux  and  the  bottom  is  sustainer 
current.  The  first  three  msec  show  a  dc  output  where¬ 
upon  the  flux  begins  to  oscillate.  At  this  point  there 
is  a  floating  zero  level  due  to  detector  sensitivity. 
The  lower  figures  are  expanded  time  scales  at  incipient 
breakup  of  flux. 
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Figure  3.  Mach  Zender  Interferograms  of  the  Flow.  The  output  flux 
trace  is  shown  as  the  lower  curve  in  the  upper  right  with 
positive  downward.  The  flux  is  stable  up  to  about  2  msec 
whereupon  the  Q-switching  occurs.  The  three  interfero¬ 
grams  are  before,  at,  and  after  onset  of  full  scale 
oscillators.  Each  picture  shows  two  exposures  separated 
by  25  psec.  The  large  scale  waves  evident  at  4.8  msec 
are  a  result  of  the  oscillations  and  not  the  cause  as 
they  do  not  appear  earlier.  At  2  msec  notice  t:  ho 
presence  of  a  standing  transverse  wave  of  a  few  fr  inner, 
and  about  3  or  4  cycles.  It  is  this  wave  that  produces 
the  oscillations  in  the  flux. 
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In  this  section,  we  will  describe  a  model  of  flux-gain- 
acoustic  (density)  coupling  which  admits  an  instability  and 
derive  an  analytic  expression  for  the  unstable  root  in  certain 
limiting  cases.  A  subsequent  paper  will  present  a  more  rigorous 
formulation,  comparison  with  experiment,  and  experimental  methods 
of  removing  the  instability. 

B.  THE  MODEL 

A  schematic  model  of  the  C02  CWEDL  is  shown  in  Figure  4. 

The  flow,  anode-cathode,  and  optical  axis  directions  are  along 
the  x,  y,  z  axes,  as  shown  in  Figure  4. 

The  anode  is  below  and  the  cathode  above,  the  plane  of  the 
paper.  Furthermore,  we  take  the  plane  z  =  0  as  lying  midway  be¬ 
tween  the  mirrors.  Assuming  that  the  change  of  the  medium  is 
negligible  in  a  round  trip  photon  time  (2  L/c)  we  may  write  the 
following  integral  equation  for  the  complex  amplitude 


u  (x,  y,  z  =  0,  tQ  +  — )  = 

[Kop  {x'  x';  y'  y’r  2  =  °/  s/  tQ)  u  (x',  y',  z  =  0,  tQ)  dx’dy' 

(13) 

In  Eq .  (13)  u  represents  the  amplitude  for  either  the  left 

or  right  going  flux;  the  propagation  operator  Kop  is  a  function 
of  the  state  of  the  medium;  and  s  denotes  the  path  taken  by  a  ray 
from  (x',y',z'=0)to(x,  y,  z=0). 

In  the  presence  of  small  perturbations,  we  make  the  decompo¬ 
sition 


K 


op 


=  K 


op  (o) 


+  K 


op  (1) 


(t) 


for  the  propagator  and 


(14) 


u  =  exp  {i(kQz  -  wot)}[uQ(x,  y,  z)  +  ^(x,  y,  z,  t)  )  (15) 

for  the  complex  amplitude,  where  kD,  are  the  wavenumber  and 
angular  frequency  of  the  light.  both  uQ  and  u,  are  assumed  to 
vary  slowly  over  a  distance  comparable  to  a  wavelength.  By  defi¬ 
nition  of  u0 ,  K0p  and  of  what  constitutes  a  laser  mode  we  have 
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u  (x,  y,  z  =  0,  t  +  — )  = 
o  1  o  c 


'op  (O)  (x'  x<'  y'  y' '  2  =  0)  UQ  (x* ,  y',  z*  =  0)  dx’dy’ 


Keeping  terms  to  first  order  in  Eq.  (13)  and  using  Eqs.  (14)  -  (16) 
we  have 


[»o<'o  +  v1  +  U1  (to  +  V>]  exp  <-  “'o' c'>  1  ’J' 

+  /K op  (II  uo  dlt'a>"  *  JK op  (o)  ' 


K  ,  ,  u  dx'dy* 
op  (o)  o  1 


u.  dx'dy’  (17) 
op  (o)  1 


If  2Lu)0/c  =  niT  where  n  is  an  integer  (i.e.,  we  have  a  longitudinal 
mode)  then 


Ul(x,  y,  2  =  0,  tQ  + 


=  JK l(to)uo  dx'dy'  +  jKo  ^lt0)dx'dy 


Expanding  the  left-hand  side  of  Eq.  (18)  in  a  Taylor  series,  we 
obtain 


Ul(x,  y,  z  =  0,  t)  +  — 


ux  «  J  Kx  uQ  dx'dy '  +  J  Kq  ux  dx’dy' 


2L  _ 
c  3 


[K  -  6(x-x')  6  (y-y')l  u.  (x',  y',  z'  =  0,  t  )dx'dy' 


rr  “  j  IKo  -  «<*- 
+  /Ki  ,Xl  x': 


y ,  y  * )  u  dx’dy' 


In  physical  terms  the  second  term  on  the.  right-hand  side  of 
Eq.  (29)  represents  scattering  "into"  Uj,  whereas  the  first  term 
is  leakage  "out  of"  uj .  If  u^  is  an  approximate  eigenmode  of  the 
loaded  resonator  with  eigenvalue  close  to  unity,  the  first  term 
on  the  right-hand  side  will  be  second  order  and  can  be  dropped. 

We  investigate  the  meaning  of  this  condition  in  greater  detail  in 
Appendix  A.  We  have  then 
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j Kx (x,  x ' ,  y. 


y  '  )  u  dx  1  dy ' 
o 


In  most  applications  it  is  the  flux  perturbation  which  is 
experimentally  of  interest  and  \vhich  couples  to  gain  and  density 
perturbations  in  the  medium.  To  obtain  the  flux  equation  corre¬ 
sponding  to  Eq .  (21)  we  multiply  both  sides  by  uc*  and  take  the 

real  part  of  both  sides.  We  obtain 


It  (uo*ul  +  uoUlM  =  ~k 


We  note  that 


[uo*  j  K ^  uq  dx'dy'  +  c.c.  (22) 


=  +  *1  =  luo  +  U1>2  =  UQ2  +  (uo*ui  +  VV*  +  0  <  1  U1  I  2 ) 


so  that 


$1  =  u0*ul  +  u0ui*  (24) 

We  now  assume  that  the  lowest  mode  uQ  is  independent  of 
x,y  (i.e.,  spatially  uniform),  so  that  uQ  may  be  taken  outside 
the  integral  on  the  right  hand  side  of  Eq.  (22).  This  will  be 
most  valid  for  small  equivalent  Fresnel  numbers,  but  is  a 
reasonable  first  approximation  in  any  case. 


o  [  /Ki  dx'dy '  +  c,c’] 


9*  _  CV1 

at  at 


As  we  have  noted  will  be  directly  related  to  the  perturbations 
of  the  medium.  We  now  turn  to  expressing  K ^  in  terms  of  density 
and  gain  perturbations. 

Assuming  the  medium  reacts  slowly  in  a  time  2L/c,  the 
propagator  KQp  can  be  written  in  the  form 

Kop  (x,  y;  x*,  y';  z  =  0,  s,  t)  = 


rx ,  y 

exp  | i  ~  n(x,  y;  x',  y ' ;  t,  z  '  =  0,  s)  d  sj 

iY'  fCj  0t 
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Where  the  path  is  the  line  segmented  geometric  path  from  (x',  y ' ) 
to  (x,y)  and  n  is  the  complex  refractive  index.  An  illustration 
for  the  typical  points  is  given  below  in  Figure  5. 

For  the  unperturbed  medium,  we  have,  simply 

/•U,  y) 


op  (o)  2 XL 


/•ix, yj 

Xl  exp  {  1  X  /  no  ds  } 


(x' ,y' ) 


Hence,  if 


n  =  n  +  n. 
o  1 


"op  (1) 


K  ,  .  =  K  , 
op  (o)  op  ( 


o)  H^x/hi  ds}  -1] 


For  sufficiently  small  values  of  n^  (less  than  one  fringe  dis¬ 
placement  on  an  interf erogram)  we  may  expand  the  exponential  in 
Eq.  (29)  to  obtain 


K  ~  ,•  2_rr 

op  (1)  X 


r  p(x,y) 

I  n1  ds 

L  Jt  X,  I  ,,  I  \ 


(x' ,y ' ) 


'op  (o) 


Substituting  Eq .  (30)  into  Eq.  (25)  there  results  for  either  the 

left  or  right  running  flux 


9<|)1  _  C(^o  2iri 

at  2L  x 


•  (x,y) 


'op  (o) 


Nr'-"1  '  l  f 

n^  ds  |  dx ' dy ' 
(x ’ y  ' )  ' 


'  fKop*  (o)  if 


<x,y) 


n^*  ds|  dx'dy' 
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Figure  5.  Illustrating  the  Propagator  via  which  the  Medium  at 

(x1 ,y ' )  on  the  Nth  Pass  Influences  the  Light  Amplitude 
at  (x,y)  on  the  (N  +  l)st  Pass 
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where  Bg  is  the  Glndstone-Dale  constant  for  the  medium  and  .  a 
is  the  density  of  the  gas  mixture  at  one  atmosphere  pressure. 
Noting  that,  by  definition 


n 


R 


and 


we  have,  finally 


(37) 


(38) 


x';  y,  y ' )  Pj  (x',  y',  t)  dx'dy5 

y  /  y’)  g-|_  (x‘,  y',  t)  dx'dy '1  (39) 

To  appreciate  the  physical  content  of  Eg.  (39),  we  transfor 
spatially  as  follows: 

(a)  In  the  y  (anode-cathode)  direction  we  take  a  finite 
Fourier  cosine  transform. 


3<J>1  C(+>o 

tt  <*■  y-  u  -  ir 


f 


aI(x,  x’; 


(b)  In  the  x  direction  we  perform  an  exponential  Fourier 
transform.  Since  the  steady-state  flux,  t0,  is  nonzero  only  over 
a  finite  region,  we  get  convolutions.  To  finesse  this  problem 
we  smooth  in  kx~space  over  a  region  2/b,  where  b  is  the  length 
of  the  flow  channel.  That  is,  we  define 


.K  +l/b 
x 


IK  2/b 


<f>1  (K^)  dx 


'K  -1/b 
x 


We  may  now  cast  Eq.  (39)  in  the  form 


(40) 


d*lm,k 

dt 


c*o  X"' 

2L  L—4  (aRjnm',KK’  °lm'K'  +  aimm'KK'  glm'K,)  (41) 

m  1  ,  K  1 
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We  now  postulate  that  one  acoustic  mode  will  be  favored  above  all 
others,  and  that  this  mode  has  no  variation  on  a  scale  loss  than 
b  (in  the  flow  direction)  and  has  m  =  l,  say,  in  the  y  direction. 
The  reason  for  this  assumption,  as  well  as  a  prescription  for 
estimating  £,  will  be  given  below. 

Equation  (41)  now  becomes  the  pair  of  equations. 


1  Cf 

0,0  _  _ o 

2L 


dt 


£  ,  o 


i  R  o 

(_  00,00 


+  2;m  gl 


0,0. 


(42) 


(where  is  the  medium  thickness) 
M 


d<f 


1 0  c 

£  ,  o  _  o 


1 


£ ,  o 


R££,oo  Po 


dt  2L 

In  Eqs.  (42)  and  (43)  we  have  tacitly  assumed 


+  2£m  91 


£ ,  o 


(43) 


a 


Imm1 ; oo 


mm 1 


(44) 


since  the  main  effect  of  a  gain  variation  is  to  change  the  flux 
in  the  same  spatial  manner.  We  note  further  that  Eq.  (42)  is  an 
equation  for  the  rate  of  change  of  the  output  flux  and  as  such 


R 


0)1,00  1,0 


(45) 


can  be  considered  to  be  a  time  varying  modulation  of  the  cavity 
coupling.  Equation  (40),  on  the  other  hand,  describes  the  varia¬ 
tion  of  the  "£"th  Fourier  component  of  the  flux  distribution. 

The  following  physical  interpretation  can  be  given  to  Eq .  (40) . 

When  the  laser  is  turned  on  there  will  be  some  acoustic  noise 
(density  fluctuations)  between  anode  and  cathode  which  can  be 
expanded  in  a  Fourier  series.  Each  acoustic  mode  constitutes 
a  phase  diffraction  grating  for  the  incident  flux  iQ.  In  Appendix 
B  we  will  show  that  the  following  statements  can  be  made  about 
the  scattered  light: 

a.  The  acoustic  gain  md  optical  modes  will  have  a  trans¬ 
verse  variation  with  period  a/n,  where  "a"  is  the  anode-cathode 
separation,  and  n  is  given  by 


n2  =  a2/AL  (46) 

zi  * 
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b.  A  piano  wave  scattered  by  an  acoustic  wave  will,  after 
one  trip,  produce  a  flux  oscillation  with  the  same  transverse 
wavelength  as  the  original  acoustic  wave. 

c.  The  variation  of  the  acoustic  perturbation  in  the  optical 
propagation  direction  can  be  neglected. 

Furthermore,  we  note  that  if  we  take  uQ  as  a  constant  for 
the  lowest  mode,  then  the  higher  Fourier  components  u^;  (recall 
‘PlU  =  u0*  ult  +  uo  ul?,*)  are  orthogonal  to  uQ  and  it  is  reasonable 
to  assume  that  the  higher-order  modes  are  representable  by  sines 
and  consines  to  the  same  extent  that  the  lowest  order  mode  is 
representable  by  a  constant.  This  simple  picture,  incidentally, 
allows  us  to  make  a  realistic  estimate  of  the  influence  coefficient 
aR U,oo  (henceforth  referred  to  as  a)  in  Eq .  (43). 

Technically,  the  quantity  a  enters  into  the  model  equations 
as  the  influence  coefficient  from  a  sinusoidal  density  variation  on 
one  pass  to  a  sinusoidal  flux  variation  on  the  next  pass  of  the 
laser  radiation. 

The  radiation  diffracts  through  the  standing  acoustic  wave, 
which  acts  as  a  phase  grating,  and  completes  a  round  trip,  rein¬ 
forcing  the  density  wave.  Thus  reinforcement  arises  by  having 
the  perturbed  radiation  in  the  same  physical  shape  (same  spatial 
Fourier  component)  as  the  density  wave.  The  reinforcement  is 
accomplished  via  heating. 

The  perturbation  in  the  amplitude  of  the  scattered  light 
(due  to  diffraction)  is 


u 


1 


1)  u  ~  i  u  6 
o  o 


(47) 


Where  uQ  is  the  incident  amplitude  and  6  is  the  (sinusoidally, 
spatially  varying)  change  in  the  optical  path.  The  change  in  the 
flux,  <{>  r  is  given  by 


VI  *u.  +  u  u  *  ~ 
o  1  o  1 


(48) 


Hence,  in  one  pass 


(49) 
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Whore  Hm  is  the  optical  length  through  the  medium,  Fg  is  the 
Gladstone-Dale  constant,  and  pa  is  the  density  of  the  medium 
at  one  atmosphere 


The  equation  in  which  a  enters  is  of  the  form 


c  Ml£ 

-TT1  “  II  <2  Vn  +  “  <501 

From  this  equation,  the  change  in  one  pass  due  to  a  density  per¬ 
turbation  is  given  by 


2L/C  TT 


By  comparing  Eqs.  (49)  and  (51)  we  obtain 

2ttR  p  a 
~  q  o  m 


We  now  proceed  to  a  consideration  of  the  gain  and  heating 
equations  necessary  to  complete  the  model.  The  gain  equation  for 
the  flowing  C02  medium  is  (fs  22  saturation  flux) 


/la-  +  v 

yat  3x  J 


_  2  _  -2<L  4  p 
T  T  d)  _ 


Where  t  is  the  effective  deactivation  time  from  the  upper  lasing 
level  of  the  CO2  molecule,  including  the  effects  of  the  nitrogen 
in  the  gas  mix,  and  P  is  the  pumping. 

Once  again,  we  consider  the  perturbed  form  of  Eq.  (53) 

(v  is  the  mean  flow  velocity) 


4  v  -- 


V’o 


Vi 


The  effect  of  the  flowing  term  is  '  lamp  the  growth  of  any  per¬ 

turbations  with  a  damping  time  of  t.iu  order  of  the  convection 
time  through  the  cavity.  This  is  of  Hie  order  of  b/v  where  b  is 
the  beam  dimension  in  the  direction  of  flow.  Hence,  Eq.  (54) 
becomes 
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(55) 


3gi  _ 

rgi 

Vgl 

go(r-l) 

Where  r 

4> 

5 1 *♦* 
s 

at  ~ 

T 

b 

T  *0 

The  £th 

component 

of  Eq . 

(55) 

will 

be  of  the  form 

8gl<>  +  V 

rgl £  go(r-1}  *11 

at  b  gl£ 

x  x  <tQ 

(56) 


To  complete  the  coupled  set  of  equations  we  use  the  fact  that 
perturbations  in  the  gain  and  flux  amount  to  a  perturbation  in 
the  heating  profile  in  the  medium.  This,  in  turn,  will  affect 
the  density  perturabations ,  completing  the  loop. 

The  acoustic  equations  for  a  flowing  medium  and  non- 
isentropic  flow  can  be  obtained  from  th<'-  equations  of  momentum, 
continuity,  and  from  the  first  law  of  thermodynamics.  In 
Appendix  C  we  show  that  the  heating  equation  is  of  the  form 


Where  Pj  is  the  perturbation  in  the  power  per  unit  volume  given 

by 


P1  =  T  (go*l  +  gl^o)  + 


gl*s 


(58) 


where 


v  =  v-,  -  v 


1' 


(59) 


and  V3,  vj  represent  the  upper  and  lower  frequencies  of  the  lasing 
transition.  Combining  Eqs.  (57)  -  (59)  with  the  prescription  (Cf. 
Appendix  C) , 


9  v 

..  'v  „ 

3x  b 


(60) 


(where  b  is  the  dimension  of  the  beam  in  the  direction  of  the 
flow)  we  obtain 
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)[(£♦£)  -  cs2  v2]<’i  " 

lY-i)  V2  j(-i  (g0+!  +  g^0)  +  —  9^3 


For  the  S-th  component  we  have 


2  2 
-  3 

y 


Where  "a"  is  the  anode-cathode  spacing.  Hence, 

+  [(it  +e)2  -  cs2  (t)2]“u 

=  {^L)  (t)  if  (go*U  +  9U*o>  +  if  9U« 


(go*u  +  gu*o>  +  —  gu*s 


Equations  (43),  (56),  and  (63)  constitute  one  set  of  coupled 

equations.  To  find  the  normal  inodes  of  the  system  we  use  Laplace 
transforms  of  our  variables  p^,  and  <J>^ 


?U(S)  -  ♦^(t) 


e  dt 


etc.  This  allows  us  to  make  the  replacement 


j£  ♦1(t)  -*■  S  <)>1(S)  ,  etc. 


In  matrix  form  our  coupled  equations  become 


h^o 

A  g1£/gc  -  o 


pU/po 


JtT  * 
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Where 


-  eg  -r- 
^ o  L 


r  v 

S  +  -  + 
x  b 


i  -»(¥)2  ^ 


(S+& 


/  \  2  2  /  \  2 

(s+b)  +cs(V) 


Equation  (66)  contains  both  the  normal  modes  and  the  instability 
criteria  for  the  system.  The  normal  modes  can  be  found  from  the 
secular  equation 


A  (S)  |  =0 


Equation  (68)  can  be  rewritten  in  the  form 


Sn  =  0 


where  yq  arc  functions  of  the  system  parameters.  The  roots  of 
Eqs.  (68)  and  (69)  are  the  normal  modes.  An  instability  obtains 
if  one  of  the  roots  has  a  real  part  greater  than  zero.  The  con¬ 
dition  that  no  root  has  a  real  part  greater  than  zero  can  be 
formulated  as  a  condition  on  the  coefficients  yn  (the  Routh- 
Hurwitz  conditions).  We  now  proceed  to  a  determination  of  the 
stability  boundary  (Eq.  (98)). 

C.  STABILITY  BOUNDARY 

Equations  (67)  -  (69)  yield  the  following  values  for  yn : 


'  Sv'SA2  +  SRA2||Sb2  *  Srb2>  -  |S  (Y  -  1> 


(v)2 
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The  Routh-Hurwitz  conditions  require  that  for  stability  the 
following  conditions  obtain: 


Yn  >  0  (n  =  1...5) 


Y4Y3  ~  Y2Y5 


Yo  -  Y . 


V  4Y3  "  Y2Y5/ 


Y0Y5 


/  Y2Y5  ] 

\  Y°V3  -  -77-; 

)  /  Y  4 Y 1  -  Y 0Y 5\ 

Y2  •  ’ 


To  make  the  algebra  tractable  we  take  the  leading  terns  in 
£qs.  (70)  -  (75).  In  determining  the  leading  terms,  we  note  that 
for  typical  systems 

SKB2  »  SRA2  »  Sv2'  SB2  IB5> 


Then  we  take 


v  e  c  2  2  P  q 
Y0  ~  v  bRA  bRD  t 


7  7  7  7  v  1 

Y,  »  Snn  Snn  +  2S  S„n  +  Pa  — 
1  1  RA  RB  V  RB  V 


Y2  *  3SV  W 


Y3  *  SRB 


.  .'AVISO  [  vtnrTi 


where 


/ ,  \  2  q  1 

p  s  - » (.;)  ^ 


Now  Eq.  (84)  is  equivalent  to  the  condition 


2  2  2  2  3 

Y0Y2Y3Y4  +  2  *’oYlY4  ‘  YftY,  Y  /,  +  Y,  YoY„  ( Y-,Y  ~  "  Y,  Y-  > 


0 ' 3  '4 


1 ' 2 ' 4  '  3  1  4 


1  '  4 


Inserting  Eqs.  *86)  -  (91)  in  Eq .  (93)  and  using  the  con¬ 

ditions  of  Eq.  (85)  we  obtain,  after  some  tediou  s  a  1  (j  t '  L>  j”  <Li  t 


{  o  r  2Pot  S  S  v. 

'JSV  +  2Sb»  4SV  sb  sra  srb  +  - T— L-  (1  +  3  v  V> 


> 


Or,  if  a  <  0, 


Pa  < 


2<V>  SRA2  SKB2 


1 


1  +  3(SVD  ~ 


The  Eqs.  (81)  -  (83)  arc  trivially  satisfied  except  for  the  co 
dition  that  >0.  If  a  >  0,  this  impos.es  the  condition  that 


Pu  <  (Si) 


RB 


2 


( 


Combining  Eqs.  (95)  and  (96),  we  have 


-2 

P 


S  i 
v 


]  +  3S  i 
v 


V 


< 


a 


P  i 
✓  -Jii- 


C-  2  q  2 

- RA  ‘  RB_ 
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as  our  stability  criterion. 
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Utilizing  the  definitions  contained  in  Eqs.  (76) 
and  Eq.  (92),  we  have,  alternatively 


-  2 


2L  _v_ 
b  *  C„ 


where 


’A 


1  +  3  —  • 

v 


2u0  P 

-  <±  — ^  -2  £  (y-1)  <  . 

.  X  p  m  4  c 

vi  a  S 


-  (80) 


2Tj  v 
b  *  C 


S 

(98) 


H  Po  V 


(99) 


D.  DISCUSSION 

Equation  (98)  can  be  put  in  a  more  useful  form.  We  note 

that 


^S  ^SA  Ip 


(100) 


T  = 


*A  — 
po 


(101) 


^A  =  0 


o 


AA  p. 


(102) 


Where  t^,  f‘SA,  denote  the  values  of  i,  ‘g,  ; ^  at  a  pressure 

of  one  atmosphere.  For  a  value  of  iR  <  0  (for  which  t  he  left  side 
of  the  inequality  of  Eq .  (98)  obtains)  the  limiting  value  of  .  /. 

which  we  denote  by  (p0/p/\)  Cr  will  be  given  by  °  “ 


(103) 
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Let  vis  consider  a  typical  system  with  the  following 
paramters 


P 


A 


0.6  kcj  m 


v/Cs  =  0.5 


^SA 

X 


2 

.  ,  .  ,  ,,  ,  2 
1.6  x  I  0  V* / m 

10  ^  m 


( 'i  -  1 )  =  0.5 

6  =  1.7x10“ 

g 

Cg  =  360  m/ sec 

b  =  7  cm 

.  5 

t  =10  sec 

o 


(105 


Under  these  conditions  the  maximum  allowable  density 

(pressure)  is  given  by  Eq.  (103)  as  ( ,-■  /:■)  -  0.38. 

o  A  c  r 

It  is  interesting  to  note  that  Eq.  (98)  does  not  contain 
the  i lux  explicitly  and  is  thus  independent  of  the  level  of 
pumping.  We  must  bo  careful ,  however,  to  note.'  1 5v.it  E  ; .  (98) 

holds  only  in  Eq.  (85)  holds.  This  implies  a  condition  on  t  )v-  • 
flux  exceeding  some  minimum  amount.  Also,  Eq.  (^5)  as sumes  1 5 tv' 
substitution  of  Eqs.  (86)  -  (91)  for  Kqs.  (70)  -  (75)  .  If  the 

latter  are  used,  ;Q  will  enter  the  stability  condition  explicit! 

The  linear  instability,  if  it  exists,  will  as  we  have  se.-i 
involve  an  exponent i a  1 1 y  orowinn  aeourtic  wave.  This  wive  nee  i 

h.ive  a".  1  v  . t  v  :  ■  to  r  -  -•  v.  1  t  in  : 

We  me.  v  see  t  :  s  oo.nl  it. it  *.  y  :  v.'s  is.  . : .  o  nr.  i  :  . 

will  t  ere.  to  ho  oleo  ■  t  -  •  '  •  ( any  •-  vt  :  it  -  ns::.-  '  •'  ' 

they  have  a  small  can:  in;  time  com.!  a:  <  v  to  tie-  lain-  ’  i  u  x  i  •  to 
(i.e.,  Sv  <<  SB) .  Hence, 


1  f 
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and 
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<t> 


o 


a 

2L 


_1_ 

£  7T 


[21  g,  „ 
l  iim1£ 


(107) 


and  the  amplification  can  be  quite  significant  (in  practice 
~  103  -  104)  . 


Equation  (98)  suggests  difficulties  in  trying  to  run  a 
high-powered  CO2  EDL  in  a  cw  mode.  A  solution  suggested  by  this 
analysis  is  to  damp  the  acoustic  mode  further  by  means  of  an 
acoustic  absorber:  for  example,  an  acoustic  muffler  on  the  cathode. 

This  has  been  tried  and  found  to  be  successful . (2 , 3 )  jn  a 
subsequent  paper  we  will  investigate  the  effects  of  an  acoustic 
absorber,  as  well  as  more  realistic  kinetics  and  compare  our 
results  with  experiments  undertaken  at  AERL. 


2.  Kellcn,  P.F.,  Mattson,  A.C.,  Ahouso,  D.R.,  and  Voder,  M.J., 

Optical  Engineering  340  (1979). 

3.  Korff,  D.,  Glickler,  S.I,.,  and  Daugherty,  J.D.,  "Acoustic 
Instability  Model  for  High  Power  cw  EDI.  Lasers,:  30th  Annual 
Gaseous  Electronics  Conference,  Palo  Alto,  Cal.,  Oct.  18-2-1, 

1977  . 

jr  * 


.  ' AVCO  1  V!  Ml  n 


SECTION  III 


ALIGNMENT  SENSITIVITY 


In  the  previous  section  we  have  ianored  a.iiii  variations  in 
the  direction  of  flow.  Such  variations  can  lead  to  extreme  align¬ 
ment  sensitivities,  as  we  shall  now  show. 


There  arc  two  effects  involved.  The  first  i  s  the  variation 
of  output  coupling  with  position  of  the  optical  axis  on  the  feed¬ 
back  mirror.  This  is  basically  a  diffraction  effect  (4)  and  is 
illustrated  in  Figure  6. 


In  Figure  6  NCq  is  the  system  equivalent  Fresnel  number, 
defined  by 


N  =  M-r_! 

eq  2  XL  M 


(108) 


where  a  is  the  collimating  mirror  diameter,  M  is  the  system 
magnification,  and  the  other  symbols  have  bee n  explained  above, 
c  is  defined  as  the  normalized  position  (one  dimensional)  of  the 
optical  axis  on  the  feedback  mirror  (>•  =  1  ,  0,  -  1  correspond 

to  top,  middle  and  bottom  of  the  mirror,  respectively). 

The  second  effect  is  the  variation  of  the  system  gain  with 
position  of  the  optical  axis.  This  effect  arises  because  of  the 
variation  of  the  gain  with  temperature  for  a  particular  value  of 
J,  the  initial  angular  momentum  of  the  radiating  molecule.  This 
variation  arises  through  the  partition  function.  In  fact, 


ga  in 


2J  4  1  -BJ  (J  4  l )  /t 


(109) 


where  B  is  the  vibrational  constant  for  the  transition  and  T  is 
the  temperature. 

Our  flux  equation  (for  the  variation  of  the  total  flux, 
not  a  Fourier  component)  is  (x  is  the  flow  direct  ion,  and  x  =  0 
is  the  optical  axis) 


34’  j 

Tt 


C  V 


(g 


al 


V 


(110) 


4.  Kc  lion,  P .  ancl  Fmith,  M.  ,  Opt.  Eng .  ,  1 8 ,  pp.  IbV-lGO 

(March/April  1979). 
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where  u i ,  proportional  to  the  tilt  angle,  is  defined  by 


U 


1 


? 

3x  P1 


(x) 


0 


=  dimensionless  coupling  constant 


(111  ) 


We  shall  derive  below. 

The  gain  equation  along  the  optical  axis  (or  the  average 
system  gain  across  the  beam,  which  can  be  shown  to  the  equivalent) 
satisfies  the  equation 


dJl  =  _  2  _  (Vl  +  qlV  +  _2_g_  liil 

9t  t  t  <f>  dy,  dt 

s  1 


(112) 


In  Eq.  (112)  the  last  term  on  the  right-hand  side  repre¬ 
sents  the  effect  of  the  gain  variation  with  optical  axis  position. 
We  rewrite  Eq .  (112)  in  the  form 


3g.  (g  <t'i  +  lit1  )  c0bg  ?;i, 

^1  cj  ^ovl  ^To  2  Jo  _J__1_ 

at  t  1 4>s  po  at 


(113) 


where  a2  is  a  dimensionless  constant  whose  value  will  be  found 
below. 


The  heating  equation  takes  the  form 
3t  Ot2  -  Cs2  ax2)  P:(x)  =  (Y  -  1 )  a2f  (x)  (go<J1  +  g^o)  (114) 


Inhere  f  (x)  denotes  the  variation  in  gain  along  the  flow  direction 
(gQ  and  g^  are  assumed  to  have  the  same  spatial  dependence).  The 
equation  for  ui  (h  pj '  (0))  then  becomes 


-  C 


a .  )  p,  ’  (x) !. 


=  (y  -  i)  dvJftx) 


x  =  0 


(g. 


A 

h  l 


c-?l 


where  the  prime  denotes  derivatives  with  respect  to  x. 
Henceforth,  we  take 


(115) 


(116) 
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to  simplify  the  analysis.  If  a  standing  wave  exists  in  th>-  flow 
channel  (analogous  to  an  open-ended  organ  pipe)  then  L'.q.  (llf>) 
follows  naturally.  In  any  event,  it  serves  as  a  reasonable  ap¬ 
proximation.  Similarly  we  take 


(y  -  1)  3x  f(x)|x  =  Q 


<-  (Y  -  1) 


(117) 


where  03  is  derived  below. 

Our  three  coupled  equation  for  this  case  are  given  by 

/*1  (s)  \ 

A  (  g"i  (s>  )  =  ° 

(s)  ' 


(118) 


where  (r  =  1  +  7—  and  v  is  the  gas  flow  velocity) 


-  c<J> 


v  r 

s  +  t-  +  — 
b  t 


a  jb 

eg  <J>  - 
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(119) 


(Y  "  1) 


3  Yo 


4>  (Y  -  1)  "  (s  +  f)  (s  +  £)*  + 


A  Routh-Hurwi tz  stability  criterion,  using  leading-order  terms 
as  in  the  previous  chapter,  leads  to  the  stability  requirement 


where 


.  _  ,  2  1 

(y  "  ala3  ~  r  ’  ~  C- g~ 

Ys  c  s^o 


_  3  .  b 

-  p  C  and  =  - 
rA  o  s  c  v 


(120) 


(121  ) 


We  shall  now  derive  expressions  for  and  (which,  as 

we  shall  see,  is  ~  «2> . 
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To  determine  ^  we  must  relate  changes  in  c  to  changes  in 
p'  (x  =  0) .  We  have 


2JL  *  .  ii  „  i  , 

A  9tilt  A  ^g  p  ix  M 
J  a 


(122) 


where  ? >,j  is  the  thickness  of  the  optical  medium  and  Itilt  is  the  tilt 
of  a  ray  due  to  refraction.  The  optical  axis  will  then  tilt  by 


A*  =  °tiit 


(123) 


Furthermore , 


~2  he  R^Aijj 


(124) 


where  d^  ,  R-^  are  the  diameter  and  radius  of  curvature  of  the 
feedback  mirror,  respectively.  Hence, 


4e  =  M  .  6  .  ,  .  erm 

d,  M-l  Rg  M  p. 


(125) 


Now,  let  us  denote  the  change  in  gc  when  c  changes  by  (1/Neq)  as 
ng0-  Then, 

Ag  =  Ac 


=  (ngQ)  Neq  Ac 

=  tngJ  •  - 


o  eq  d.  M  M  -  1  p. 


Ap'  (0) 


(126) 


alyob 


Ap' (0) 


Therefore , 


X,  _JM _  R]_  J4  ^o 

al  1  Neq  M-l  d1  l'g  Pa 


For  a  typical  case 


^  ~  n  i  m  -  1  20  =  200 

n  «  o.i,  Ncq  =  io,  m  =  2,  —  --  -0--T 


(127) 


=  0.1, 


2  x.  10~4,  -----  =  -  -  50 

d  ^  2  cm 
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u  < 


so  that 


«  a 

We  note,  furthermore  that  nj  will  be  much  less  than  its 
alculated  value  if  the  system  is  aligned  at,  e.g.,  points  A  or 
'  in  Figure  6.  (The  calculation  in  Eqs.  (126)  and  (127)  is 
for  point  B  in  that  figure.) 

We  may  schematically  denote  the  variation  of  system  gain 
with  optical  axis  by  the  diagram  in  Figure  7. 

An  analysis  identical  to  that  above  leads  to 


n  .  ..  4M  R  ^o 

a2  Neq  M  -  1  b  d.  g  p 

1  ci 


(128) 


Now,  ot2  arises  from  the  temperature  variation  of  the  gain 
as  given  in  Eq .  (109), 

Taking 


T 

inlet 


-  17  3°K , 


=  380° 


we  obtain 


a2  ~  0.5 


for  a  typical  case. 

Since  02,  013  arc  both  dimensionless  constants  due  to  the 
same  effect  we  take  them  to  be  equal.  We  then  obtain  the  follow¬ 
ing  stability  boundary. 


,  .  .2  f  4MVm _ l2  ,  2  /^o\ 

n  (Neq)  [(M  -  l)bdj  bg  \  pa  J 
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1  c 


For  a  numerical  example,  tak 
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cq 

2  0m 

b  = 

1  0  cm 

0.05 

eM  “ 

1  m 

-  4 

2  cm 
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T  = 

C 

10  P  sec 

SYSTEM  GAIN 


In  such  a  case,  we  have 

>  4  (unstable) 

2b  3 

An  instability  in  this  system  would  result  in  large  wander¬ 
ings  or  the  pointing  angle  and  possibly  to  line  jumping  if  the 
optical  axis  moved  sufficiently  far  so  that  the  gain  were  greater 
for  a  different  J  number. 

Even  in  the  absence  of  an  instability,  of  course,  the 
matrix  A  determines  the  system  response  to  external  perturbations 
In  particular,  if  we  have  mirror  vibrations  ( ac  perturbation)  the 
coupling  ensures  that  we  would  obtain  a  response  at  the  acoustic 
frequencies.  Such  a  response  would  not  obtain  in  standard  treat¬ 
ments  . 


Finally,  we  stress  that  if  the  system  is  aligned  at  a 
stable  point  (A  or  A'  in  Figure  6)  no  instability  obtains.  >1 
(or  n)  is  negligable  and  Eq.  (129)  is  a utomat ica 1 ly  satisfied. 
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SECTION  IV 


MODE  MEDIUM  INTERACTIONS:  DIFFERENTIAL  FORMULATION 


A.  INTRODUCTION 

The  analysis  in  this  section  will  deal  mainly  with  an  attempt 
to  recast  the  instability  discussed  in  Section  II  in  terms  of  more 
conventional  medium  instabilities  such  as  are  found  in  plasma 
physics.  Some  of  these,  such  as  the  free  electron  laser,  harmonic 
generation,  and  optical  modulation  by  sound  waves  are  beneficial. 
Others,  like  thermal  blooming,  se 1 f- focusing  and  mode-medium 
interaction  (MM I )  are  deleterious. 

The  analysis  in  the  previous  two  sections  emphasized  the 
Huygens-Fresnel-Kirchhof f  integral  solution  to  the  wave  equation. 

In  this  section,  instead,  the  differential  equations  will  be 
analyzed  and  comparison  with  the  other  derivation  will  be  made 
along  the  way.  The  differential  formulation  has  the  advantage 
that  it  can  be  easily  extended  into  the  nonlinear  regime,  in  which 
the  amplitude  of  the  perturbation  has  built  up  to  a  large  value. 

It  also  has  the  advantage  that  it  demonstrates  how  the  instability 
discussed  in  Section  II  (hereafter  referred  to  as  MMT-1)  is  related 
to  more  conventional  instabilities.  To  do  this,  Maxwell's  equations 
including  a  time-varying  conductivity  (to  account  for  cavity  and 
medium  losses)  and  a  time-varying  dielectric  function  will  be  used. 

It  will  be  seen  that  MMI-1  results  when  the  conductivity  is  modulated 
and  when  all  the  electromagnetic  waves  have  the  same  unperturbed 
frequency.  This  is  a  form  of  nonresonant  amplitude  modulation  and 
instability,  which  is  different  from  the  conventional  method  of 
amplitude  modulation  in  which  resonant  coupling  of  electromagnetic 
modes  which  are  separted  by  the  modulation  frequency  occurs.  On 
the  other  hand,  Brillouin  scattering,  which  is  primarily  a  phase 
modulation,  results  when  the  dielectric  function  is  modulated. 

The  coupling  is  strongest  when  the  electromagnetic  modes  are 
separated  by  the  modulation  frequency.  This  is  the  mechanism  that 
will  bo  studied  in  this  section.  The  equations  of  motion  for  the 
electromagnetic  and  acoustic  waves  are  derived.  The  secular 
determinant,  is  solved .  It  is  shown  to  have  unstable  root  s  in  the 
limit  in  which  the  ratio  of  cavity  flux  to  the  saturation  flux  is 
small,  so  that  gain-flux  oscillations  can  be  neglected.  Finally, 
the  nonlinear  regime  is  briefly  examined. 

B.  MODE-MEDIUM  INTERACTIONS 

Our  starting  point  will  be  Maxwell  '  s  egu.it  ions  with  th<‘ 
appropriate  boundary  conditions  (cavity  losses)  and  mod i nr  g  li n 

9 
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and  losses.  Since  the  laser  is  in  an  unstable-resonator  config¬ 
uration,  the  boundary  conditions  are  rather  complicated  (see 
Appendix  D) .  Thus  to  simplify  matters,  we  will  lump  all  the  losses 
in  a  phenomenological  conductivity  function.  The  medium  gain  will 
be  computed  semiclassically  from  the  imaginary  part  of  the  suscep¬ 
tibility.  A  simple  model  of  this  is  discussed  in  Appendix  E. 

Thus,  the  governing  equations  are,  in  mks  units, 


V  •  D  =  0 


(130) 


B  =  0 


(131) 


V  x  E  =  - 


(132) 


V  x  H  =  J  + 


(133) 


where  D  =  eE,  J  =  oE,  c  contains  the  effects  of  dielectric  polari¬ 
zation  and  o  contains  both  the  medium  losses  and  the  cavity  loss. 
From  Eqs .  (130)  -  (133)  the  following  wave  equation  can  be  derived 

for  transverse  electromagnetic  waves: 


V  E  -  y 


°dt2 


(134) 


(Uq  is  the  free-space  permeability) .  Assume  that 

e  =  c  (0)  (r)  +  e  (1)  (r,t) 
o  =  o(0)  (r)  +  (r  ,t) 


(135) 

(136) 


The  entire  scattering  calculation  can  actually  be  done  in  terms  of 
a  or  t: .  The  only  reason  for  keeping  both  is  to  stick  to  the  con¬ 
vention  which  associates  the  gain  of  the  system  with  the  imaginary 
part  of  c  and  the  losses  with  the  real  part  of  c.  Similarly, 
scattering  is  associated  with  the  real  part  of  <■  and  not  the 
imaginary  part  of  a  ev^n  though  formally  (with  the  proper  inter¬ 
pretation)  it  can  be  calculated  in  terms  of  either  function.  In 
what  follows  we  maintain  the  convention.  Substituting  Eqs.  (135)  - 
(137)  into  Eq .  (134)  leads  to, 


( 0  )  ,  .  3  2  E 
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3t. 


(0)  TE  _ 
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0  L. 
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This  equation  will  now  be  solved  in  a  perturbative  fash  i J 
However  rather  than  expanding  in  normal  modes,  we  will  use  the 
multiple-time-scale  expansion  of  secular  cquat.  i  < -ns  .  ((■ )  Kxia.udi: 
E  =  EO  +  E(l)  +  ?A2K  realizing  that  F.(0)  :  0  since  there  is  no 

steady-state  space-charge  field,  leads  to 


V2E(1> 


(0)  o2e  _  (0)  ?k{1) 

U0C  2  U0°  3t  ~  0 

at 


(138) 


V2E(2) 
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(139) 

If  e  ^  and  in  turn  depend  on  the  electric  fields,  then  third, 

and  higher-order  equations  will  result.  Equation  (138)  describes 
the  equilibrium  eigenmodos  of  a  loaded  resonator  including  all  the 
losses.  This  equation  for  the  unstable-resonator  coomet ry  is 
discussed  in  Appendix  D  and  in  more  detail  in  another  put  of  this 
report.  The  normal  modes  of  E  ( 2 )  are  the  same  as  those  for  i;U) 
(since  the  left-hand  sides  of  the  equations  are  identical). 

However,  since  E  (2)  is  driven  by  products  of  f  i  rr.t  -order  guard  it  ic 
it  is  conceivable  that  the  left-hand  side  is  driven  at  a  natural 
frequency.  If  this  should  happen  to  be  the  case  then  our  pertur¬ 
bation  expansion  will  break  down.  The  reason  is  that  ];(••)  contain 
a  secular  response  which  will  make  E  ( 2 )  >  e(1).  To  get  around,  tin 
problem  those  components  of  the  right-hand  side  of  E ( 2 )  are  incor¬ 
porated  into  the  equation  for  F,U)  ,  which  is  now  allowed  to  have 
a  slowly  varying  amplitude.  First  the  e i uonnodes  of  pq.  (lib)  liav 
to  be  found.  Towards  this  end  let, 


E(1)  (r , t)  = 


(1)  -io)kt+ikz 

e  +  c  •  e 


(140) 


where  c-c  stands  for  complex  conjugate  and  P 
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where 


,22  (0 
,  k  e  .  q 

1  2  aR0N0  +  1  k 


(142) 


g  ,  gc  are  the  equilibrium  gain  and  loss  terms.  The  loss 
includes  medium  and  output-coupling  losses;  it  is  related  to 
o(0)  by,  gc(0)  =  (r)ok/k.  The  relationship  between  the 

dielectric  function  and  the  medium  susceptibility  is  given  in 
Appendix  E.  The  dipole  moment  per  molecule,  g i von  by 
p  =  (aE  +  c-c),  results  in  a  macroscopic  polarisation  P  =  UqP- 
Finally,  since  d(D  =  c(0)e(1)  =  ipEd)  +  F ,  we  find  that  in  oi 
simple  semi-classical  picture  c(0)  =  e  q  ( 1  +  t  0  -  ig(0)/kj. 
The  equation  for  the  equilibrium  gain  is  given  by 


V  .  g^V0)  = 


4  SI - 


(0) , (0) 


(143) 


where  v'  is  the  equilibrium  flow  velocity  of  the  gas,  ■  is  the 
collisional  deactivation  time  of  the  upper  state  and  P  is  the 
pumpinq  rate.  The  solution  to  Fqs.  (142)  and  (143)  constitutes 
the  equi  1  ibri.um  eigenmodes  of  the  resonator,  since  the  resonator 
boundary  is  included  in  c  (0)  .  The  details  are  given  in  Ap: ondix 
D.  Any  perturbation  in  the  equilibrium  fields  can  be  expanded  in 
the  complete  orthonormal  set  of  eigenfunctions  (which  are  assum'd 
to  exist).  However,  since  a  secular  response  is  ant  i  o  i  i  ate  : ,  the 
Fourier  coefficients  will  be  assumed  to  bo  time  dependent  and  the 
appropri a te  set  of  equations  describing  the  time  evolution  of  the 
coefficients  will  now  be  derived.  That  is,  let 


E  ( r ,  t )  =  X  uk  ( t )  Fj,  (r)e 


-  ii.  ^t  +  i  k  z 


(144) 


and  sill's,  t  i  tut  e  this  into  Fq.  (137),  keeping  in  mind  th.it  Ur  (t)  is 
slowly  varying  and  Ek  ( 1  )  satisfies  Eg.  (1*11).  The  pstml  at  ion 
expansion  in  Kqs.  (138)  and  (139)  thus  is  redone  to  take  into 
account  anticipated  secularity  in  Die  response  of  k(2).  Assume 
also  that 


r!l)(r,t)  -  V  (t,  (n(r)tri'hl,lk:'- 
k 
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0(1)  ( r  ,  t )  -  X  h.  (t)  -  (1)  (r)c 

k  k  k 


-!,kt4lkZ 


(3  4G) 


where  >k  (r)  and  'p.  (r)  are  the  appro;  riate  spatial,  dielectrii 

and  loss-modulation  e  i  qon  functions  to  be  determined  se  1  f-cc.  :  i  s- 
tontly.  Note  that  in  F.qs.  (144)  -  (140)  k  is  a  summation  ird  ; 
there  is  no  implication  that  the  k's  are  th.e  same  in.  these  equa¬ 
tions  .  Substituting  Egs.  (  144)  -  (  146)  into  I’q .  (137)  yields, 
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Notr'  +hat  the  secular!  ty-produci  ng  terms  lead  to  slow  variations 
in  the  mode  amplitudes  as  expected.  To  pick  out  a  particular 
eigenmode  on  the  left,  multiply  the  left  and  right  sides  by  the 
complex  conjugate  of  the  left  and  integrate  over  the  resonator 
volume  and  time: 


3uk  ic2  ^  f 

g-jm-  =  +  2 “  X  /  dvdt  exp  [i  (wk-u>k  , -u>k„  )  t-i  (k-k ’-k"  )  z] 

k  k  '  ,  k  "  J 


x  (wk,  +  o;k„)  l1oak,uk"Ek  ^  Ek"  ^ck' 


+  -  X  /dVdt  exp  [  i  (o,  -w,  , -w,  „  )  t- i  (k-k  ' -k  "  )  7.  . 
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2  ( 1  )  *  ( 1 ) 

where  the  eigenmodes  are  normalized  such  that:  /d  ri:^  J'k 
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There  are  several  possible  cases  of  interest' : 

(i)  Consider  first  the  case  in  which  there  is  no  dielectric 
modulation  d)  (r)  -  0)  . 

Taking  the  limit  that  (.:k  ,  (the  slow  modulation  of  the 
gain  and  loss  term)  compared  to  optical  frequencies  is  negligible 
requires  that  .g.  =  w^** .  That  is,  the  two  elect  roman  net  i  c  waves 
have  the  same  natural  (unperturbed)  frequency  of  oscillation. 

One  way  this  could  happen  would  bo  if  k"  ~  k,  i.e. ,  the  lasing 
mode  is  "coupled"  to  itself  by  loss  modulation  arising  from  the 
acoustic  wave.  Alternatively,  the  two  modes  may  be  distinct, 
although  degenerate  in  frequency;  this  case  is  treated  in  detail 
in  Section  II.  In  the  self-coupling  ease  we  obtain 


3u}.  2 

3~=c;i0(bk.e  )u. 


"fd2r  i: 


.  (1)*,  (l)r  (1)  n.q. 

k  LV  Lk"  (149) 


with  a  similar  equation  for  Uj.,, .  Making  the  following  identifica¬ 


tion 


aRfp*0 


(r)  rk  ,  (1 '  (r)  i:k„  (r) 


results  in  the  equation  used  in  MMI-1  with  ;e  being  the  density 
perturbation  by  the  acoustic  wave.  The  gain  term  can  bo  obtained 
from  the  polarization  term  as  is  done  in  Appendix  D.  The  above 
identification  of  nj,.  >  (1)  with  density  fluctuations  is  not  unrea¬ 
sonable  since  the  cavity  and  medium  losses  for  simple  systems  can 
be  related  to  the  quality  factor  of  the  system  which  in  turn 
depends  on  the  index  of  refraction.  Besides  the  case  vy.  =  ..k„  it 
is  also  possible  to  consider  the  resonant  triplet  ■  •  =  ’.q: "  +  <■  p  *  , 

k  -  k'  +  k"  where  now  the  two  electromagnetic  modes  are  separated 
by  the  acoustic  frequency.  The  only  difference  is  that  now  the 
-i^k .  t 

e  term  in  Eq.  (149)  is  replaced  by  unity.  The  right-hand 

side  is  again  related  to  the  density  fluctuation.  This  is  the 
case  of  conventional  amplitude  modulation  by  a  resonant  triplet. 

(ii)  Consider  now  separately  the  case  of  no  cavity  and  medium 
loss  variai  ion  (  '  v,  *  0)  but  no*  >  include  a  dielectric  modulation 

f p •  ¥  0.  To  proceed  further  wo  develop  a  simple  mode  1  for  •  y » . 
The  specific  modulation  of  the  medium  that  is  to  be  studied  is 
that  due  to  an  acoustic  wave. 


7.  M.  Born,  E.  Wolf,  Principles  of  Optics,  Pcrgnr.vm  Press,  MY  (199^) 
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From  microscopic  considerations  (see  Appendix  E)  it 
follows  that  the  dipole  moment  of  a  molecule  is  given  by 
P  -  ("I’eff  +  C'C)  where  a  is  a  frequency-dependent  proportionality 
constant  and  K0ff  =  F.  for  a  qas.  The  macroscopic  polar  i  zat  ion  ( 8  ) 
of  the  medium  is  then  given  by  (for  N  molccules/cubic  centimeter) 

ff  =  Np  (150) 

with  the  dipole  moment  given  by  p  =  (rtE  +  c-c).  Note  that  a,  U 
both  depend  upon  the  intensity  of  the  field.  Since  the  secular 
part  of  the  response  comes  from  products  of  first-order  quantities, 

USin<!  “sec  “  cI1)e<1)  =  t0('“ORNk(1>  -  1  TT  '  '’h'1’  +  C" C  leads 
to 


8_uk 
3t  . 


+  cUk„/d2rgk(1>(t,r)Ek*(1»Ek„'1> 


(151) 


where  it  was  assumed  that  the  fast  acoustic  exponential  variation 
has  been  extracted  via  the  requirement  that 

w  =  mk,  +  «k„  (152) 

k  =  k'  +  k”  (153) 

and  the  fast  transverse  variation  in  and  g^d)  is  taken  care 

of  by  the  electric  field  variations  which  have  different  transverse 
mode  structure 

v  =  ki  +  V  (154) 

where  k  j_’  is  fixed  by  the  acoustic-resonator  geometry.  In  the 
plane  wave  case  the  spatial  integrals  can  actually  be  done  and 
lead  to 


3u 


k 
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6q  =  ft  q  r  /  m  is  the  Gladstone-Da  1  c  coefficient,  i  a  is  the  atmo¬ 
spheric  density  and  *  denotes  complex  conjugate  quantit.i  es. 
Equations  for  cj y  t  will  be  derived  in  the  next  section. 

The  only  term  which  can  drive  the  gain  equation  resonant  1' 
[that  is,  so  that  Eqs.  (152)  -  (1  54  )  are  satisfied]  is  the  flux- 
term.  Thus,  the  secular  part  of  the  gain  equation  is 


Vk- 


^0  T 


(r-1)  ukuk" 


(157) 


where  r  =  1  +  Co/\J's  and  -  0  -  ccquv^/2  and  (-,  a  =  3/-.  +  V0/b  (whom 
V0  is  the  flow  velocity  and  b,  the  anode-cathode  distance)  in,  th* 
total  damping  of  the  gain  term  including  convection  and  deactiva¬ 
tion  of  the  upper  level. 

The  equations  describing  the  medium  are  the  continuity, 
momentum  anci  heat-transpert  equations, 


-^r  +  V  •  pV  =  0 


(158) 


|r  +  v  •  w  -  -  -'X 

3 1  p 
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|“  +  V  •  VP  =  f-te  t  V  .  Vp)  +  ( Y 


(160) 


where  p,  V,  P,  dQ/dt  are  the  mass  density,  fluid  velocity,  pressure 
and  heat  deposition  rate;  y  is  the  ratio  of  specific  heats. 


Linearizing  Eqs.  (1  58)  -  (160)  would  produce  tin. 
acoustic-wave  equation.  However,  since  the  heat  deposit' 
can  have  components  at  the  acoustic  frequency,  this  lead: 
secular  response  in  the  acoustic  equation.  This  secular 
once  again  taken  care  of  by  letting  the  am.jilit.uder;  be  sic 
varying.  The  linearized  acoust  i  c-wave  equation  (for  the1 
density)  ir;  given  by 
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where  Vq  is  the  transverse  flow  velocity  and  Ca2  =  y RT  is  the 
acoustic  velocity.  The  heating  torn  has  been  included  in  antici¬ 
pation  of  the  secular  i  ty .  The  heating  term  is  obtained  for  a  D'.i 
laser  by  conservation-of-energy  arguments;  it  takes  the  following 
form ; 


(162) 


where  vj/vo  is  the  ratio  first  excited  state  to  ground  state  and 
C'C  stands  for  complex  conjugate.  Since  P  does  indeed  vary  at 
the  acoustic  frequency  (note  that  use  of  Eqs .  (152)  -  (154)  was 
made)  it  will  produce  a  secular  response  in  Eq .  (161).  The  slowly 

varying  piece  of  the  acoustic  wave  is  produced  by, 
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where  y  a  is  the  convective  damping  of  the  acoustic  waves,  ■  is 

given  by  Eq.  (162)  and  pj^  >  =  mNj^<  where  m  is  the  mass. 


The  linear  instability  is  obtained  by  assuming  that  one 
of  the  electromagnetic  waves  has  large  amplitude  and  hence  decays 
into  the  second  electromagnetic  wave  and  the  acoustic  oscillations. 
Thus  assuming  the  vq  '■  constant  and  using  Eqs.  (  1 56)  ,  (157)  and 

(163)  leads  to  the  determinantal  equation  for  the  instability. 
Assuming  that  u^„  -  eot,  ,  ~  e^*t,  g^  ,  ~  e  -'-**1  leads  to, 


Setting  the  determinant  of  Eg.  (164)  equal  to  r.-ro,  leads  to  the 
polynomial  equation  for  the  instability.  Defining 
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where  SQ  =  Gladstone-Dale  constant  and  p  is  atmospheric  density, 
we  obtain  a 


(6)  (<5  +  Yg)  (6  +  Ya)  +  a1(6  +  ya)  "  i«2(6  +  >  )  +  i8  =  0  (165) 

The  acoustic  damping  is  given  by  ya  =  VQ/b  where  VQ  is  the  trans¬ 
verse  gas  flow  velocity  and  b  is  the  anode-cathode  distance.  The 
damping  in  the  gain  equation  is  given  by  ',-g  =  VQ/b  +  1  / x  where  t 
contains  all  the  collisions!  relaxation  rates. 


C.  NONLINEAR  SATURATION 


If  now,  it  is  assumed  that  u^  does  not  stay  constant  but 
decreases  significantly  (as  energy  is  transferred  to  the  other 
waves)  ,  it  is  necessary  to  include  an  evolution  equation  for  uj. . 
That  is  Eq .  (155)  must  be  included  together  with  Eqs.  (156) , 

(157),  and  (163).  Unlike  the  three-wave  equations  of  plasma 
physics,  which  can  be  solved  in  terms  of  Jacobi  elliptic  functions, 
this  set  of  equations  can  only  bo  solved  numerically.  Some  insight 
into  these  solutions  can  be  obtained  by  examining  the  conserved 
quantities.  Using  Eqs.  (155)  -  (157)  and  (163)  leads  to, 
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This  can  further  be  reduced  by  using  the  equation  for  g,  ,  . 
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The  convective  damping  has  been  included  simply  by  replacing 
V01/~X  in  the  acoustic  equation  by  VQ/b  where  b  is  the  anode- 
cathode  spacing. 

In  order  to  understand  Eq .  (167)  it  is  necessary  to  examine 

Figure  8.  During  the  linear-lasing  stage  the  second  electromag¬ 
netic  wave  is  negligible  and  the  acoustic  noise  has  not  made  a 
transition  in  the  anode-cathode  direction.  Once  the  primary 
electromagnetic  wave,  u^  has  boon  saturated  and  reaches  steady 
state  and  the  acoustic  noise  has  had  time  to  propagate  across  the 
anode-cathode  separation,  the  linear  regime  of  the  instability 
sets  in,  whereby  the  primary  electromagnetic  wave  feeds  the 
secondary  electromagnetic  (em)  wave  up”  and  the  acoustic  noise. 
Eventually  the  primary  cm  wave  is  totally  depleted.  This  results 
in  the  secondary  em  wave  also  dying  out  in  one  or  two  round  trips 
in  the  resonator  because  of  the  output  coupling.  However  Eq .  (167) 

implies  that  the  nonlinear  process  of  supplying  energy  back  to  the 
primary  wave  will  not  take  place  because  of  the  dissipative  nature 
of  the  acoustic  waves.  The  energy  can  be  made  to  flow  into  the 
secondary  em  wave  and  the  acoustic  wave  but  will  not  flow  back  to 
the  pump.  Thus,  once  the  pump  is  depleted  the  instability  stops. 
The  entire  process  begins  again,  however,  because  uo  builds  up 
due  to  the  laser  properties  of  the  medium.  The  acoustic  wave  may 
remain  into  this  region  because  of  its  relatively  long  decay  time. 
This  feature  may  lead  to  a  degradation  of  beam  quality  of  pulsed 
lasers,  as  well. 

D.  SUMMARY 

The  physical  processes  involved  in  the  instability  are 
illustrated  in  Figure  9. 

Consider  a  large-amplitude  wave  at  frequency  >,o  and  wave 
number  kg.  In  the  presence  of  a  density  fluctuation  at  frequency 
w  and  wave  number  k,  a  second  electromagnetic  wave  ( q  -  , 

kg  -  k)  is  excited  through  the  nonlinear  polarization.  This  wave 
is  a  higher-order  mode  of  the  unstable  resonator  and  lienee  could 
have  a  higher  loss  rate  than  the  fundamental  mode.  Therefore,  it 
may  not  be  excited  in  the  absence  of  the  density  fluctuation.  The 
two  electromagnet ic  waves  in  turn  amplify  the  density  fluctuation 
through  the  heating  term.  In  addition  they  excite  fluctuations 
in  the  equilibrium  gain  through  the  flux-saturation  term.  This 
completes  the  feedback  "loop."  Since  the  heat  deposition  if;  a 
dissipative  process,  this  instability  is  rmneonsorvat ivo  and  turns 
itself  off.  The  reason  for  this  is  that  once  all  the  energy  from 
the  primary  wave  ( te q  ,  kn)  has  been  deposi  t  ed  into  t  he  secondary 
electromagnetic  wave  and  the  acoustic  wave,  there  is  no  way  of 
returning  the  energy  to  the  primary  wave  as  would  be  the  case  of 
ordinary  conservative  nonlinear  instabilit  ies.  Pecause  of  the 
additional  loss  of  the  secondary  wave  (higher-order  mode)  it,  too, 
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Figure  9.  Physical  Processes  Involved  in  a  Three-Wave  Instability 
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leaks  out  once  the  primary  wave  has  decayed.  The  whole  process 
restarts  because  the  primary  wave  is  once  aqain  built  up  by  the 
lasing  medium.  The  acoustic  waves  may  remain  present  throughout 
this  phase  because  of  their  relatively  lony  damping  time  constant. 

The  equations  we  have  derived  in  our  differential  formulation 
reduce  to  those  of  Section  II  when  the  two  electromagnetic  modes 
are  degenerate  in  frequency,  as  they  must. 


The  process  of  dielectr 
showrn  to  be  unstable.  This 
including  the  effects  of  ga 


ic  modulation  was  studied  and  also 
is  a  form  of  Brillouin  scattering 
in-flux  oscillations. 
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SECTION  V 


UNSTABLE  RESONATORS  WITH  SATURABLE  GAIN 


A.  INTRODUCTION 

Unstable  resonators  have  found  wi despread  use  in  hi  di- 
gain  lasers.  As  early  as  1965  )  it  was  recognized  t  h  it  t  hen. 
resonators  are  desirable  for  a  number  of  reasons.  These 
include : 


1)  The  large  mode  volume  that  can  be  obtained. 

2)  The  uniformity  of  illumination  (and  hence  energy 
extraction)  seen  by  the  gain  medium. 

3)  The  relative  freedom  (compared  to  stable  resonators) 
from  undesired  higher-order  transverse  modes. 


In  the  realm  of  geometric  optics,  the  operation  of  an 
unstable  resonator  such  as  the  one  shown  in  Figure  10  is 
obvious:  The  output  is  a  plane  wave.  However,  when  diffrac¬ 

tive  effects  are  included,  modi"  properties  are  not  obvious  and 
the  calculation  is  complicated  in  the  extreme.  The  understand¬ 
ing  of  unstable  resonator  properties  that  has  evolved  over  the 
past  15  years  has  been  derived  from  the  gradual  development  of 
ever  more  sophisticated  calculation  t  echn  i  ones  .  ^  n_  * r’ )  Prior 
to  this  report,  the  study  of  unstable  resonators  had  advanced 


9.  Siegman,  A.E.,  Proc.  IEEE  53,  u.  277  (1965). 

10.  Fox,  A.,  and  Li,  T.  ,  Bell  Syst.  .  Tech.  J.  4_0 ,  pp.  453-488. 

11.  Streiffcr,  W. ,  IEEE  J.  Quant.  Elect.  QE-4 ,  pp.  229-230 
(April  1968). 

12.  Bergstein,  L.,  Appl  .  Opt  .  1_,  p.  495  (1968). 

13.  Chen,  L.,  and  Fo]sen,  h. ,  TELE  J.  Quant.  Elect.  QE-9, 
pp.  1102-1113  (November  1973). 

14.  Horwitz,  ?.,  J.  Onl  .  Lei.  Am.  6  3,  op.  1  528-1  543  (December 
1  973)  . 

15.  Moore,  G.,  and  McCarthy,  R. ,  J.  Opt,  hoc.  Am.  67,  pp.  228- 
241  (February  1977) . 
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rc  10.  Telescopic  Unstable  Resonator,  Showing 
Colli mat eJ  Output  Beam 
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to  tho  point  where  the  effects  of  diffraction  and  deterministic 
intra-cavity  gain  distributions  were  understood  .( 1  r> )  In  this 
report  we  carry  the  level  of  understandinn  one  step  further, 
and  include  the  self-consistency  problem  associated  with  gain 
saturation  as  found  in  actual  lasers. 

B.  NOMENCLATURE  AND  STATEMENT  OF  PROBLEM 

In  this  section,  we  express  the  coupled  diffraction  and 
saturated  gain  problem  in  terms  of  the  one-dimensional  strip 
confocal  unstable  resonator,  as  shown  in  Figure  11.  Comparable 
analysis  can  be  made  using  two-dimensional  rectangular  symmetry 
or  circular  symmetry  resonators,  but  we  wish  to  avoid  these 
additional  mathematical  and  computational  complexities  here. 

It  should  be  noted  that,  for  the  case  of  empty  (no  gain)  reso¬ 
nators,  the  diffractive  properties  differ  only  in  detail  among 
the  several  geometries. 

In  the  absence  of  either  a  saturable  or  non-uniform  gain, 
the  optical  properties  of  the  strip  resonator  are  described  by 
the  well-known  integral  equation  (16) 


Xf  (x) 


/  2  i  M2 

"  V  M2  -  1 


eff 


+ 


+  e 


c 


exp 


cf  f 


.  f(x’)  dx* 


(168) 


In  this  equation 


f (x)  =  the  complex  optical  amplitude  in  the  plane  of  the 
feedback  mirror 


A 

M 


F 


eff 


the  eigenvalue  corresponding  to  fix) 
geometric  magnification 
the  so-called  effective  Fresnel  number 
(a2 / A ' L) ( (M  -  l) /2) 


A'  =  wavelength  of  light-  (not  to  hi  confused  with  '  ) 


16.  Weiner,  M.  ,  Appl  .  Opt.  18,  pn.  1  8  2  1  H  1  -1  (Pune  1  ri  7  P )  .  The 

basic  resonator  cow  it  ion  is.  d« -set  i ■■  d  in  tuim- a  '  "is  ;  aper:  , 
but  this  work  is  t  he  only  on<>  (n  ext  t  ess  i  t  explicitly  in 
terms  of  the  effective  Fresno)  number. 
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2a  =  feedback  mirror  size 

l  -  the  normalized  optical  axis  offset  or  tilt.  r 
takes  on  values  from  0  to  1 . 

Equation  (168)  is  just  a  formulation  of  the  requirement 
that  the  optica]  amplitude  distribution  reproduce  itself,  to 
within  the  constant  after  making  a  round  trip  back  and  forth 
in  the  resonator.  For  example,  in  the  limit  Feff  diffrac¬ 

tion  effects  vanish  and  Eq.  (168)  takes  on  the  form 

Xf(x)  =  f (x/M)  (169) 


which  has  as  a  solution  f  (x)  =  constant  and  \  =  1. 
readily  identified  with  the  geometric  optics  result 
wave  scaled  in  size  by  a  factor  M  every  round  trip, 
output  coupling  seen  to  be 


This  is 
of  a  plane 
with  the 


L  =  1  -  (1/M)  (170) 

which  is  just  the  fraction  of  the  mode  not  blocked  by  the  feed¬ 
back  mirror. 

In  many  cases  of  practical  interest,  however  (e.n.,  infra¬ 
red  lasers),  a  geometric  description  is  not  sufficient  and  dif¬ 
fraction  effects  must  bo  included,  even  though  the  optical 
properties  may  be  largely  geometric  in  character.  This  corre¬ 
sponds  to  solutions  of  the  integral  equation  for  which 


1 


eff 


100 


1 . 5  <  M  <  4 


For  a  given  combination  of  M,  F0ff  and  f  ,  then*  are  several 
solutions,  each  with  its  eigenvalue  and  ei uenf unct i on ,  corre¬ 
sponding  to  the  several  transverse  modes  of  t  lie  resonator. 
With  the  eigenvalue  equation  defined  as  in  Eq.  (168) ,  the 
(one-dimensional)  output  coupling  I,c  ( i  )  associated  with  each 
mode  is  given  by 


Lc(i)  »  1 


(171  ) 


We  may  as 
output  coupling 


ess  the  effect,  of  diffraction  by 
or  equivalently,  i  >  1  ,  as  F(1<  (  i 


men i t  ■ 1  r i ng 
vu  r  i  ■  -d  . 
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Figure  3,  fro:n  Ref.  13,  shows  such  a  curve  for  M  =-  2.9,  ..  -  0, 
and  Fef f  between  3  and  20.  An  overall  periodicity  in  F,.ff  is 
apparent,  as  well  as  degenerate  crossing  points  (o.q.  ,  Ft.f  f  —  ft.  9) 
where  two  modes  share  the  same  value  of  |  >■  [  .  In  an  actual  laser, 
such  points  may  cause  undesirable  mode  hopping.  A  very  similar 
pattern  of  diffraction-induced  output  coupling  variation  exists 
if  Fcf f  is  kept  constant,  but  the  position  of  the  optical  axis 
(c )  is  varied . ( , 1 8 )  Even  if  the  mode'  are  not  actually  degen¬ 
erate,  the  magnitude  of  the  output  coupling  ripple  is  important 
for  assessing  mode-medium  interaction  of  foots .  This  latter  situ¬ 
ation  is  perhaps  the  more  likely  in  an  operating  device,  as  it 
can  readily  lie  induced  by,  say,  a  mirror  vibration,  or  by  a 
density  perturbation,  leading  to  f-lMI  instabilities. 

In  any  realistic  high-efficiency  laser  the  gain  will  always 
be  saturated.  The  effect  of  the  presence  of  gain  saturation  or. 
theoretical  curves  such  as  depicted  in  Figure  12  has.  hereto:  or-  • 
not  been  assessed,  and  provides  the  purpose  for  this  study.  in 
the  main  theoretical  analysis,  we  consider  only  ali  resona¬ 

tors  (c  =0),  but  in  Section  E  we  show  how  tilted  resonator 
behavior  can  be  related  to  the  aligned  case. 

C .  METHODOLOGY 

For  the  present  analysis  we  have  developed  a1.  : 

combining  a  number  of  techniques  which  have  not  1  .  ••  us-  .!  '■•••••:•  : 

before.  The  result  is  an  iterative  procedure  which  •  d  •  -s 

not  only  the  lowest  loss  mode,  but  also  1  hr'  h :  sher-o :  r  .*•••.  .  :.a- 

tor  modes,  simultaneously.  Even  for  a  highly  saturated  n>  hixm, 
convergence  is  obtained  in  five  iterations  or  less  ,  ar.d  t.her.  •  i  c 
no  difficulty  obtaining  a  solution  in  the  presence  of  node  a  r 
eracy . 


We  outline  here  the  techniques  which  compr  i  re  !:;•••  compo¬ 
nents  of  the  algorithm.  However,  since  tic  t  hoorsl  :  r  1 1  as:  >.  c 
of  these  arc  in  places  quite  intricate,  the  reader  i  s  referred 
to  the  original  references  for  complete!  detail:-.. 

1 .  Initial _ Empty  Resonator  Calculation 

For  a  given  set  of  resonator  parameters  (M,  P.  <  <  ,  with 
assumed  zero  for  the  rest  of  the  discussion)  ,  t  •  emp*  y  rswrul  sr 
eigenvalues  >  £  and  corresponding  modi'  distributions  :  ]  ( x )  ire 
determined .  The  asymptotic  theory  of  resonators  is  is  .•  !.  (11) 


17.  Perkins,  .7.,  and  Cason,  C.,  Appl  .  Phys_.  _  T,<  t  t  .  11  ,  jp.  1(18- 

200  {August  1977)  . 

18.  Kelli-n,  P.,  and  Pm  i  t  h ,  M.,  Opt. .  Fnu.  1ft,  pp.  1  r  >  7  —  1  GO 

(March/April  1  979).  '  . 
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This  theory  can  be  summarized ,  in  fairly  physical  terms, (15) 
as  follows.  The  amplitude  function  fj_(x)  is  ex 


series  of  slowly  modulated  cylinder  waves,  plum 
'Ey  as  follows: 


januee  as  a. 
a  core  term 


fi(x)  = 


N 

E 


i  k  rn(x) 


i  k  rn(-x) 


f  ■  (x) 
n ,  l 


^n(x-) 


+  f  • (-x) 

-  n ,  i 


rr  (->:) 
n 


/\ 

+  fi(::) 


(172) 


The  origins  of  the  cylinder  waves  arc  successively  further  removed 
virtual  images  of  the  feedback  mirror  edaes  (Figure  13) ,  corre¬ 
sponding  tc  light,  that  has  made  n  trips  around  the  resonator. 


Then  rn  is  the  distance  from  the  n^h  virtual  imago  to  the  obser- 


<x) ;  E  {i 


n  1 


feedback  mirror 
(172)  becomes 

plane.  In  the 

(±)1/2  . 

exp 

1  ,  .2 

2  lk(x  ~  xn 

1  w  ) 

w 

\  n  / 

n 

__  ^ 

n ,  i 


(x) 


(173) 


+  exp 


r  i  2i 

2  ik (x  +  xn ) 2 

fn(-x) 

Vs' 

n 

/x 

+  fi (X) 


The  origins  of  the  cylinder  waves  are  at  the  points 


2n 


w  =  d  (M  -  1  ) 
n 


xn  =  nM 


(174a) 

(174b) 


Use  of  Kq .  (173)  in  Kq.  ( 1  (>  8 )  1  cad:;  to  a  series  of  intecrnls  whio 

arc  solved  by  the  method  of  st  ionary  phase.  (I'M  Tins  is  vat  id  f  nr 


i  I  t.  s  l  o  r 


Fef f  >>  1,  but  has  been  r.i'.own  to  aive  reasonable  res 
Fc'ff  y.  1.  When  this  i  s  done,  the  modulation  functio 
the  core  term  1  •  are  found  to  bo  expressible  in  term; 
resonator  parameters .  The  final  result  is  a  complex  pulynomi 
in  a  parameter  u : 


•  n  ,  :  ■  •  ■  ■ 
< '  f  ('  t  la  ■  r 


19.  Horn  ,  M.  ,  and  W’ol  f,  K.  ,  Pr  i  ne  i  pj  es  of  o?>;  ire. ,  St  b  J.'d. 
(Pergamon,  New  York,  1  97S)  ,  pp.  7  \>2 -7  S4  . 
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Figure  13.  First  Virtual  Image  of  Feedback  Mirror  Ed.aos. 

Coherent  addition  of  contributions  from  numerous 
such  virtual  images  gives  rise  to  the  complicated 
mode?  structure  of  the  resonator. 
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1  =  - 


2k  \  2  F 


N  +  1 


•  £ 


n  =  1 


exp  2  k  i  F  ff  8 


e  f  f  1  n  /  '  n 


+  exp  l  2  r  l  F 


-1  V+l/2  In 


ef  f  ' n  I  n 


for  symmetric  modes  when  n  =  N  +  1 


(17  5) 


The  plus  and  minus  signs  denote  symmetric  and  antisymmetric 
modes.  Here 


e  = 


1  -  M 
1  +  m" 


(176) 


and  the  original  eigenvalue  X  is  found  from 


X  =  -4 


(177) 


The  lowest-loss  mode  corresponds  to  tlie  larges  value  of  |  X 
The  corresponding  eigenfunction  is  found  from  using 


f  ,  (x) 
n ,  l 


=  (viM)nd  ~  M  2n) 

(1  -  M-2) (a  +  M_nx) 


(178) 


y\  N  +  1  _  o  _  i 

fi(x)  =  [a(1  ~  M  )] 


2  iriF 


(179) 


in  Eg.  (172)  .  Since  the  expression  for  f ^  has  a  factor  1  -  u  in 
the  denominator,  and  since  |xj  =  ] / j n  |  ~  1  for  the  lowest-loss 
mode,  it  dominate;;  Eq .  (172)  for  regions  inside  the  shadow  boun¬ 

daries  (x  -  +  Ma)  and  represents  the  "geometric"  part  of  the 
solution . 
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The  number  of  virtual  imager.,  N,  which  are  needed  to  fix 
the  order  of  the  polynomial,  can  be  estimated  in  the  following 
way.  (20)  Consider  the  confocal  unstable  resonator  shown  in 
Figure  14  with  a  spherical  wavefront  emanating  from  the  common 
focal  point.  Radiation  within  the  angle 


0 


1 


D/M 


(180) 


which  makes  one  round  trip  through  the  cavity,  will  exit  magni¬ 
fied  to  the  whole  mode  size  D.  Similarly,  radiation  within  the 
small  core  angle 


0 


n 


(181) 


which  makes  n  round  trips,  is  magnified  M  times,  and  also  exits 
with  the  full  mode  size  D.  [This  core  region  is  where  (x)  in 
Eq.  (172)  arises.]  If  we  desire  the  cavity  to  bo  "diffraction 
limited"  in  the  usual  sense,  the  exiting  radiation  must  have  a 
cavity  angular  spread  of 


0 

c 


(182) 


By  requiring  the  angular  spread  of  the  core  region  to  match  that 
of  the  whole  cavity,  i.e.,  0N  =  0C,  we  have 


whence 


D 


MNf. 


X 

D 


(183) 


N 


In  (8  M  Feff) 
In  M 


(184) 


In  practice,  it  turns  out  to  be  viseful  to  add  two  or  three  to 
N  as  calculated  by  Eq.  (184),  to  guarantee  that  N  will  be  large 
enough.  A  similar  formula  used  by  Horwitz(l^)  is 


N 


In  (250  F  ff) 
In  M 


(181.) 


20.  Anan'ev,  Y.,  Sov.  J.  Quant.  Elect.  3,  pp .  615-81 7  (1575). 
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Intensity  Ronorma 1 i zati o n 


The  lov.’cst- loss  mode  intensity  distr  i  but  ion 


fn(x,d)  I  2 


•  0\  aiui  |  corresponciinct  to  'Q,  is  normalized  to 
ated  intensity,  Is.  To  do  this,  we  corin'  to  a  ueo; 
intensity,  I~,  using  the  theory  of  geometric:  strii 
with  saturable  gain.  (21)  We  then  require  that  t!i« 


I  ( )  - 

the  satur- 
ic  trie 

)  resonators 
.■  average 


intensity  over  the  feedback  mirror  be 
out  diffraction: 

+a 


the  same  with  and  with- 


u 

-a 


l(x)  dx 


(18G) 


This  is  to  a  certain  extent  heuristic  and  is  strictly  true  only 
in  the  limit  Fcff  -►  ®°,  but  is  a  resonable  requirement  for  the 
range  of  Feff‘s  of  interest. 

In  the  geometric  theory,,  we  begin  with  an  assumed  fi-'ld 
form  in  the  cavity;  namely,  a  modulated  cylindrical  wave  ema¬ 
nating  from  the  confocal  point,  and  a  modulated  plane  wave 
exiting  the  laser: 


E 


ik  ( z  +  L  +  d) 
e 


f  (z  ,  r ) 


g  ( r ,  ) 


(187) 


(Nomenclature  is  defined  in  Figure  15.  The  function  f(z,r)  here 
is  similar  to  but  not  the  same  as  the  f  function  of  Fiction  I.) 
The  boundary  conditions  are  that  the  amplitudes  vanish  on  the 
mirrors : 


f (0,r) 


1 


(L  +  d) 
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L  +  d. 


_ r_ 

L  + 


and 


f (d,r) 


1 

1/2  q 
d  ' 


e 


2  i  k  I, 


By  starting  with  the  time-dependent  wave  equation 


2  2  n2} 

c  V  E  = 


r»t 
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2  c  G 


’>V. 
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21  . 


Moore,  C.,  and  McCarthy, 
227  (February  1977). 


R.  , 


m 

70 


J  ._  Opt  .  pfic  .  /V’l . 


(188a) 


(188b) 


(1 

2  21- 


AVCD  !  VI  Mill 


I 


7\ 


-&U- 

i 


Jo  070 


Figure  15.  Nomenclature  for  Geometric  Strip  Resonator 
(after  Ref.  21) 
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and  neglecting  all  second-order  derivatives,  we  obtain,  using 
the  assumed  form  for  E,  the  rate  equations 


Using  Eqs.  (192,  (193),  and  (195)  in  Eq.  (191)  yields 


dH 

G 

o 

dz 

2h  ,  1 

(  I,  +  d 

\  -2H 

_  2 

1  a- 

e  +  | 

Vd  +  L  "  * 

)‘!  J 

f 

g 

(19 


Determining  a  solution  consists  of  finding  a.  value  of  fq  such 
that  the  boundary  conditions  (Eq.  (188)1  are  obtained  when 
Eq .  (196)  is  integrated.  Then  lq  =  >  M  f is  used  in  Eq .  (186) 

to  normalize  the  diffractive  result. 

3.  Propagation  Through  the  Gain _ Kou_ion 


The  renormalized  fQ(x)  is  propagated  t  liroudi  the  gain 
region,  using  simultaneous  forward  and  backward  prouaaat  ion. 
This  process  is  illustrated  in  Ficmre  16.  Consider  a  point 
(>:,w)  in  the  region  between  mirrors.  Here,  w  is  a  coordinate 
measured  from  right  to  left  away  from  the  feedback  mirror. 


The  "forward"  part  of  the  propagation  consists  of  the 
expanding  protion  of  f0(x)  reflected  into  the  cavity  by  the 
feedback  mirror,  and  thus  propag.it  es  f  ro:n  right  to  left  in  t.  h<  • 
figure.  In  the  absence  of  gain  we  'nave,  from  the  scalar  theory 
of  diffraction,  the’  complex  amplitude  at  a  point  (x,w): 


Ucx  (x'w) 


K ( x , x  x )  fo(Xj)  exp 


ik  2 
2d  ‘'l 


(19 


The  exponential  term  is  the  phase  curvature  result  i  nq  from  t  h.e 
curved  feedback  mirror.  in  the  Fresnel  approximation  K  is  the 
propagation  function  defined  by 


K(x, 


i 

XL. 


1/2 


exp 


i  k 
21, 


(18 


The  "backward"  port  ion  i  r;  the  nit  i  p'  extent  of  fri(x)  , 
prior  to  hitting  the  mirror,  evolving  )  re  k  ward  in  t  in. •.  This 
propagation  is  effected  using  the  compl  "x-con  j  uqate  propaqa  t  i  <  n 
f  unct i on : 


24  . 
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*  jr 

UBKW(X'W)  =  J 

—  OO 

K*  ( x  ,  x x ) 

W 

dx  , 

(199) 

Thus  it,  too,  propagates  from 

right  to 

left  . 

In  this  way, 

both 

waves  move  into  the  gain  region  in  the  same  spatial  sense. 

Locally,  the  saturated  gain  is  computed  using,  again,  the 
Rigrod  formula: 


G  ( x ,  w )  = 


o 


1  +  |Ucx(x,w) i2  +  |UBKW(x,w) | 2 


(200) 


Of  course,  light  reaching  a  particular  (x,w')  must  have  traversed 
gain  in  the  region  w  <  w '  .  To  account  for  this,  we  insert  in 
Kq.  (198)  a  line  integral  term 


Ggx (x , w, Xf )  =  exp 


G(f 


(201  ) 


where  is  a  vector  connect!  na  the  observing  point.  (x,w)  w  i  t  h 
each  point  ( x  i  ,  0 )  on  the  feedback  mirror. 

For  the  backward  wave,  gain  becomes  attenuation,  so  that 


gbkw 


(x , w , x1 ) 


=  exp 


G  ( >' ) 


(202) 


is  used  iii  Eq .  (199)  .  This  scheme  permits  syst  ein.it  ic  calcula¬ 
tion  of  G(x,w)  in  the  direction  of  increasing  w. 

4.  Modes  in  the  Presence  of  Gain 

The  analysis  of  Moore  &  McCarthy  P  ’  *  pc  asymptotic 

calculation  of  unstable  resonator  modes  •  •.•  a  in  f ur.et  ion 

can  be  specified.  In  their  work,  G(x,w)  was  a  per  t  u  1  a  tod  func¬ 
tion,  but  this  need  not  L>*  so  in  nenoral  ,  ami,  as.  r.liown  abovi • , 
wc  calculate  it  from  the  mode  pattern  itself. 

To  include  the  effects  of  gain  in  the  resonator  analysis 
the  eigenvalue  polynomial  [Kg.  (17b) )  is  modified  to  read 


lb 
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1  =  - 


2  \2  F 


ef  1. 


1  /  2 


N  *  1 

E 

n  -  1 


’^(a)  exp 


i  F 


ef  f  'ii 


,-1/2 

n 


+  r '  (-a)  exp  [2  i;  i  F  ff  p  ')  6  +  1/2 
—  n  \  e  i  i  n  /  n 

•(r-H 


(20  3) 


for  symmetric  modes  when  n  =  M  +  1 


The  now  functions,  (x)  ,  are  integrated  gain  funct  ion:;  vorv 
similar  to  Eqs .  (201)  and  (202) ,  except  that  the  path  of  inte¬ 
gration  forms  n  round  trips  through  the  cavity,  start  in  j  at  the 
mirror  edge  (x  =  a)  and  terminating  at  the  point  x.  Speeifical lv. 


FMx)  ^  exp 


p 

E 

V  1 


G ( X  x ( w 1 )  , w ' ]  dw' 


o 


(20-5) 


L 

/ 


G  ( X  2 ( w ' )  , w ' ]  dw' 


In  J  G  ( 0  ,  w  '  )  dw  '  J 
o 


where  the  integration  paths  are  alone  rays  with  c  c-i  d 

im"1"1)  ix-i-fta  -  i vf- 1 ! 


Xx  (  2  '  ) 


M 


1  -M 


(.’O'.  ,) 
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In  this  analysis,  wo  always  assume  that  the  1  owest- 1  ons  node- 
determines  the  intensity  in  the  (join  region,  so  that  I-'  func¬ 
tions  have  an  implicit  i  =  0  subscript. 

5 .  Iteration  to  Convo rooncc 

The  updated  f  (>:)  amplitude  distribution  is  propagated  into 
the  gain  region  as  in  step  3,  and  a  new  f  (:•:)  calculate  d  as  in 
step  4.  This  process  is  repeated  until  convergence  is  obtained , 
by  which  we  mean  that  the  change  in  the  lowest-loss  eigenvalue, 
X0,  is  less  than  some  values 


AAq  <6  (206) 

Since  \Q  ~  1 ,  this  is  the  same  as  having  the  fractional  change 
into  A 0  be  less  than  5 . 

For  a  GqL  of  5  (i.e.,  round  trip  unsaturated  intensity 
gain  is  c^O)  ,  a  convenience’  criterion  of  '  =  0.002  requires 
<  6  iterations;  for  *  =  0.0005,  the  number  is  eight  to  ten. 

There  appear  to  be  two  reasons  why  convergence  is  rapid  with, 
this  algorithm.  First,  the  saturated  gain  calculation  always 
beains  with  an  f0(x)  derived  from  an  on:  t:y  resonator  calcula¬ 
tion.  Thus  in  a  sense  if  treats  on  in  as  a  perturb. :t  i  on .  A 
second  reason  is  that  th-’  integrated  gain  f  unctions,  d  ’  (:•;),  are 
fairly  smooth,  because  they  integrate  out  fine;  scale  longitudinal 
(w-direet  ion)  gain  rippl>.  ,  and  do  not  chance  by  much  between 
iterations . 


It  must  be-  e-'.phan  i  1  that  the  iterative  al  oorit  hn  '..'re¬ 
sented  in  this  paper  is  •  ur.dnment  al  ly  different  from  tie  Fc  i 
itcrat.  i  on  HO)  generally  used  to  simulate  lasers  with  saturabl” 
gain.  (^5)  in  the  Fox-J.i  r.cb.eme ,  the  resonator  in  team  1  ■  •••  ;ua  t.  ion 
itself  is  iterated  to  convergence ,  starting  v;  i  t  h  eon'/'  o.;n  ured 
initial  mode  shape  such  as  a  pi  a no  wave  or  Gaussian.  In  such 
schemes,  the  only  know  led  »•  •  of  the  mode  survi  vi  :ic  f  rom  pass  to 
pass  is  contained  in  the  numerical  array  which  sg- •  i  f  i  •  -  s  f, ,(:■:). 
The  mathematical  man  i  pu  1  a  I.  i  on  nee  ded  to  propao t  •  •  f,  tx)  arouirl 
the  cavity  of  necessity  accumulate  *;  error,  bet  h  fro:;;  ruck  souse.  . 
as  roundoff  error,  and  from  approx  i  mat  ions  used  in  ♦  !,,  •  ore:  ,< 
tion  scheme.  The  necessity  for  maintaining  a  liially  iccvn.sto 
copy  of  f  (:<)  implies-  high  sampling  f  rogi.vnry  coriy.irt-t}  t  -  >  typical 
structure  scale  sires,  contained  i  :i  f(x),  with  the  l.-nult  th’.t 
computer  storage  r<  on  i  rc'ent:;  dictate  a  practical  up:  limit,  in 

Fcf f ,  usually  around  unity  and  certainly  baa;  than  10. 


25.  bz  i  !:1  as  ,  F  .  ,  and  S i com an ,  A.,  Ann!  .  Oil  .  It,  pp  .  1  7  4  -  1  ?■' <'  M 
(August  l'1','"-).  The  computer  a!  go:  iti-m  b  nvri  t  ■  i  in  thin 
paper  is  the  archetype  of  a  number  ol  n.  i  mul  at  ion  c  don  in 
use  today. 
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In  contrast,  the  method  presented  here  need  only  pass  the 
eigenvalue  >.Q  from  iteration  to  iteration.  Furthermore ,  since 
the  I’ '  (x)  functions  are  sensitive  only  to  the  largest  spatial 
variatio’  s  in  gain,  a  highly  accurate  propagation  through  the 
gain  region  is  not  essential. 

D.  RESULTS 

The  basic  purpose  of  this  calculation  is  to  assess  the 
effect  of  gain  saturation  on  well-known  diffractive  node  prop¬ 
erties.  To  this  end,  wo  consider  the  particularly  illustrative 
case  of  a  resonator  with 

M  =  2. 9 

F  _  =  8.892,  9.390,  and  9.863 
ef  f 

Referring  again  to  Figure  12  these  Fcff  values  correspond  to 
(in  Horwitz's  terminology)  a  so-called  crossing  point,  a  point 
of  maximum  mode  separation,  and  a  cusping  point,  respectively. 

To  choose  appropriate  small -signal  gains,  we  consider  the 
steady-state  oscillation  condition  for  a  strip  resonator  in  the 
Geometric  limit: 


1  2  GoL/(l  +  d/Is)] 

M  °  =  1 


GoL  =  2  (1  +  (I/Is) 1  In  M 


(207) 


(208) 


We  will  consider  the  three  regimes  of  interest: 

J/1B  “  1 

I/Is  =  1 
I/Xs  >>  1 

A  very  low  value  of  I  / 1 is  the  .name  as  the  cm  p  t  y - r o  s  o  n  a  t  o  r  cal  • 
culation  and  is  calculated  as  a  matter  of  court.*'  ‘  o»*  any  ot  h*-r 
case.  As  an  example  of  the  intermediate  case  we  tnk*-  G  T , 

In  M  -  1.065,  so  that  1/1,-,  =  1.  For  1/1,.  N>  1,  we  have  chosen 
G0L  -  5.0. 
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A  comparison  of  t  ho  pro;:*:.1  method  v:i  t_h  other  resul  f  c.  is 
frustrated  by  the  paucity  of  other  recalls,  and  by  t  .fact  tbit 
since  other  methods  ur><>  V:  i  i  t  oral  i  or.  requiring  lev.-  r..rr,  :n 

results  are  available-  at  the  conditions  listed  at  the  bouinni •: 
of  this  section. 

Nonetheless,  wo  have  ca  1  n;  1  a.  tod  a  case  for  com.par  i  son  , 

with  results  obtained  by  Ken  sell  and  Chester.  (26)  j»,  tie  ir  work, 

ct  Fox-Li  strip  re  son  at  or  calculation  was  modified  to  include  a 
saturable  gain  sheet  positioned  rinht  at  the  orirrirv  mirror: 
the  amplitude  wave  as  prop.s in t ed  to  the  primary  mirror,  ampli¬ 
fied  and  saturated  us  i  ng  the  Hi<mxi  formula  (with,  local  inf  -  n- 

sity  doubl  ed  to  account  for  incomi  ng  and  outgo i  :kj  waves )  and 
propagated  back  to  the  feedback  mirror.  The  resonator  param¬ 
eters  are: 

M  =  2.5 


F 


eff 


0.64 


G  L  -  4.583 
o 


For  this  value  of  FCf  f  time  asymptotic  method  is.  be  1  cw  its  red 
of  proven  validity.  Also,  Ron  sc  h  and  CL-cter  dc  net  rrevi  vis- 
eigen  values  but  only  intensity  distributions.  Th.us  a  comparison 
of  their  results  with  ours  must  of  necessity  be  somewhat  specu¬ 
lative  . 


Figure  17  compares  bare  resonator  intensity  patterns ,  and 
gives  a  clue  as  to  how  the  asymptotic  method  break'  down  at  very 
low  FCf  f :  the  fine  scale  detail  is  lost,  though  lh"  cress  shave 
is  basical  ly  tho  same.  The  situation  with,  sain  is  compare  1" 
(Figure  18)  .  It  will  be  noted  that  there  is  a  .client  dit  fen-nee 
in  the  location  of  the  .side  maxima.  In  our  ca  1  a:  1  a t:  i  ■.  :: ,  tiny 
occur  at  slightly  higher  values,  of  x/n.  We  sue  pent  tbit  thus- 
is  due  to  an  imprecision  on  Re  ns-,  eh  and  CL-stor's  pari  in  calcu¬ 
lating  the  effective  Fresnel  number,  s  i  nee  with  discrete  i  nl  e- 
qrution  schemes  there  is  always  an  uncertainty  in  just  where  the 
feedback  mirror  ends. 

For  the  principal  ca  1  cu  1  at  i  on  s, ,  *  L  gain  was  con  f  :  n-w;  tm 
a  region  from  z  0.23  L  to  z  0.75  ],,  sir.ee  ir:  many  h  i  ?h  ;.  ve  r 
devices  (where  urisl.abl  e  resonators  find,  i  !:•  ; :  •  :  •  •  ct  u:  •  )  t  h. 

minors  atm-  c.f'mowh.at  removed  f :  o:rt  t  lasing  m"dium.  ip  ve :  , 
this  is  not.  essential  to  the  abtorithn,  and  thr  wain  can  be 
placed  arbitrarily  between  the  mirrors. 
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Results  for  F(.f  f  =  8.092  and  G0L  :  1.065  are  shown  in 
Figures  19  through  24.  Figures  19  and  20  show  intensity  and 
phase  for  the  first  pass,  i . o. ,  before  any  gain  saturation, 
but  after  renormalization.  The  units  of  intensity  are  such 
that  Is  =  1 .  The  results  after  5  passes  appear  in  Figures  21 
through  23.  The  first  depicts  f^ix)  for  various  n.  As  n 
increases,  the  curves  approach  a  limiting  form.  As  one  would 
expect,  the  r^(x)  have  maxima  and  minima  inversely  related  to 
those  seen  in  the  intensity  plot.  Inspection  of  the  bare  and 
saturated  intensity  profiles  shows  that  the  fine?  scale  diffrac¬ 
tion  features  are  very  nearly  t he  same  in  both  cases,  save  that 
the  overall  envelope  of  the  intensity  is  made  more  uniform.  Ke 
note  that  it  would  be  hard  to  claim  that  the  diffractive  com¬ 
ponent  of  the  pattern  is  exactly  suppressed  by  gain  saturation . 
Figure  15  depicts  the  saturated  gain  in  the  resonator.  The 
optical  axis  coordinate  is  in  percentage  of  w/L;  i.e.,  0  is  the 
feedback  mirror  location  and  100  the  primary  mirror.  The  gain 
scale  is  in  units  of  gain  per  optical  axis  increment.  Thus  the 
unsaturated  gain  is  seen  to  be  1.065/50  ~  0.02. 

Figures  25  through  28  are  for  Ff,ff  =  8.892  but  with  0oI, 
increased  to  5.0.  The  most  noticeable  effect  of  the  increased 
gain  is  that  the  very  small  amount  of  "stray  liciht"  well  beyond 
t  he  geometric  mode  size  becomes  significantly  amplified.  Since 
real  devices  often  have  scraper  apertures  spooi  f  .i  cal  3  y  to  (ex¬ 
clude  such  radiation,  its  appearance  in  the  present  calculation 
is  probably  not  important.  One  can  also  discern  that  the  ranee 
of  intensity  excursions  from  a  mean  value  has  boon  still  furtlv- 
flattened.  Interestinaly  enough,  the  phase  profiles  are  nearly 
identical  among  the  three  cases. 

Figures  29  throuqh  38  depict  results  for  Fpg  -  9.39,  and 
those  for  Feff  -  9.863  are  shown  in  Figures  39  through  48.  it 
will  be  seen  that  the  same  qualitative  observations  can  be  made 
for  these  cases  as  well. 

The  preceding  plots  do  not  hint  at  the  dramatic  effect  of 
gain  saturation  on  eigenvalue,  and  hence ,  output  coupling.  In 
Figure  49  a  portion  of  the  empty  cavity  plot  of  Figure  12  if; 
reproduced  for  Feff  between  0.5  and  10.5,  alone  with  results, 
from  the  present  calculations.  It  will  im  seen  that  the  mode 
degeneracy  at  ~  8.892  is  lifted,  even  for  the  lessor  value 

of  G0T,  --  1.065  corresponding  to  I  / 1 --  1.  Calculation;;  at 
closely  spaced  intervals  near  F(.f  f  ft .  ft  g  2  hav  vri  f  i<  d  the 
the  degeneracy  is  truly  lifted,  and  not  just  moved  to  a  r.r-.ny 
Fresnel  number .  It  is  also  evident  that  even  at  it."  cusp  inn 
point  the  ripple  in  mode  separation  is  reduced ,  thouch,  sum  i- 
singly,  there  is  not  as  much  difference  letwo.-n  t  la-  cy,,  ]  .Of.', 
and  GC,I,  ■  5.0  cases  as  there  is  between  t  hese  and  the  en.pt  y- 

resonator  car.e. 
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Figure  25.  i-,'  (>:)  Profile  for  c;c)T,  -  5.0.  F,,ff  :  8 .  8‘»2 . 

Compare  1  to  F inure  21,  thi  .  curves  th> 

much  higher  nain  in  the  qoonn -t  r  i  o  shuOov.- 
region . 
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Figure  32.  Intensity  Profile  for  GqL  =  1.065.  Feff 


9.39. 


'  AVCO  FVFRFTT 


NTU9 


PASS 

FEFF 


3.000 

9.390 


GOL 

M 


1.065 

2.900 


r Avca  r:vmr-TT 


GRMMq 


o 

o 


GAMMA 


.  avco  rvt  nm 


PHRSE 


.  avco FVFnrTT 


INTENSITY 


2*40.00 


(I 


8  9  10  II 

JS068  ^eff 


Figure  49.  Effect  of  Saturable  Gain  on  Mode  Degeneracy.  The 

solid  curve  is  a  portion  of  Figure  12,  representing 
an  "empty"  resonator  in  the  absence  of  gain  satura¬ 
tion  . 
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In  a  qualitative  sense,  saturable  gain  seems  to  produce 
an  effect  not  unlike  that  of  mirrors  with  rounded  edges  (27)  or 
tapered  ref  1  activity ,( 28)  which  reduce  the  influence  of  diffrac¬ 
tion  on  mode  properties.  This  is  in  spite  of  the  fact  that,  as 
we  have  seen,  the  intensity  profiles  themselves  are  net  really 
smoothed  on  a  fine  scale. 


Horwitz  has  pointed  out  that  empty  resonator  mode  degen¬ 
eracy  occurs  for  Feff's  less  that  some  critical  value,  with 
that  critical  Feff  increasing  quite  rapidly  as  M  approaches  1. 

In  other  words,  the  propensity  for  mode  degeneracy  is  emphasized 
for  low  Fcff  and  low  M.  We  hypothesize  that  for  each  value  of 
Fc££  and  M  corresponding  to  a  mode  degeneracy,  there  is  also  a 
critical  GQL  which  will  lift  that  degeneracy.  Though  we  have 
not  explored  this  thesis  exhaustively,  we  have  performed  a  num¬ 
ber  of  further  calculations  using  G0L  =  2.  For  M  =  2.9  there 
turn  out  to  be  no  degeneracies  for  FGff  2.  When  M  is  reduced 
to  2.0,  degenerate  mode  behavior  is  suppressed  for  Fcff  •’  13. 

For  the  empty  resonator  case,  the  corresponding  critical  Feff's 
are  17  and  43,  respectively. 

This  result  is  of  considerable  importance  for  device 
design,  as  it  frees  Foff  in  many  cases  of  interest  from  being 
constrained  to  half-integer  values,  and,  by  inference,  permits 
a  more  relaxed  view  of  resonator  alignment  tolerances. 

In  light  of  the  present  results,  one  is  tempted  to  specu¬ 
late  how  the  mode  structure  of  circular  mirror  resonators  would 
be  affected  by  gain  saturation.  For  such  resonators,  it  has 
been  shown  that,  in  the  absence  of  gain,  the  in-phase  diffrac¬ 
tive  contribution  from  the  mirror  perimeter  loads  to  strong 
mode  degeneracy,  persisting  well  into  what  might  be  thought  of 
as  the  regime  of  geometric  optics. (29)  it  is  certainly  not 
intuitive  whether  gain  saturation  could  be  expected  to  lift,  such 
a  strong  degeneracy.  Clearly  this  is  an  area  which  warrants 
further  study. 


E.  IMPLICATIONS  FOR  MM I 


In  an  empty  resonator,  the  peak-to-peak  excursions  in 
the  value  of  the  \0  eigenvalue  are,  in  the  region  whore  mode 
degeneracy  exists,  about  equal  to  the  geometric  optics  limit 
symmetric  mode  separation.  (See  Figure  12  or  the  figures  of  _ 
Refs.  14  and  16.)  Since  as  Fc f f  -►  •«,  -»  1  and  11  ^  v  ~  l/M*2->. 

This  amounts  to  a  fractional  variation  in  the  output  coupling  of 


27.  Santana,  C.,  and  Felsen,  L. ,  Appl .  Opt.  17 ,  pp.  2239-2243 
(15  July  1978). 

28.  McAllister,  G.  ,  et_al._,  1  F.F.F.  J  .  Quant  .  F.  1  cc  t .  QN- 1  0  , 
pp.  346-355  (March  1974). 


29. 


Butts ,  R . ,  and  Avi 
1078  (August  1978) 


J .  Opt. . 


hoc .  Am 68  ,  pp .  1072- 
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ALC  \20  -  \\ 


o _ n  ~  2  ax 

M  -  1  ~  M  -  1 


which,  for  M  -  2.9  is  about  25%. (30)  if  saturable  gain  effects 
are  included,  on  the  other  hand  this  variation  is  clearly  cut 
at  least  in  half,  depending  on  the  value  of  G0L. 

Even  though  these  results  are  derived  for  varying  Feff 
with  e  =  0 ,  we  now  present  an  argument  that  they  apply  equally 
well  to  the  varying-e  case,  with  Fef£  constant,  because  the 
diffractive  processes  involved  are  essentially  the  same. (31) 
Consider  the  tilted  resonator  shown  in  Figure  50,  with 
e  =  gq  ^  0.  We  want  to  consider  two  effective  Fresnel  numbers, 
derived  from  the  portion  of  the  resonator  on  either  side  of  the 
optical  axis: 


(210a) 


F(2)  =  F  (0) 

*eff  eff 


(1  -  e) 


(210b) 


if  g  is  small  compared  to  unity,  then 


(211a) 


(1  -  2c) 


(211b) 


Thus  as  g  increases  from  zero,  the  two  half-resonator  Fresnel 
numbers  move  in  opposite  directions  away  from  F^-f.  In  parti¬ 
cular,  the  output  couplings  associated  with  each  half-resonator 
follow  the  usual  periodic  curves,  as  indicated  in  Figure  51,  so 
that  when  2  c  f(Q1  is,  an  integer,  f(3)  and  F^l  are  each  an 


that  when  2  c  F^j^  is 
integer  away  from  F^jr 


The  observation  that  the  next-to- lowest- 1  or 
a  geometric  optics  asymptote  of  ~  m^.25  )ia< 
previously,  but  is  apparent  from  a  close  e: 
1 X (  vs  Fcff  plots. 

A  more  rigorous  exposition  of  the  contribul 
effects  of  the  two  edges  can  be  found  in  C 
L.  Felscn,  Apnl.  Opt.  17 ,  rp.  23r'2-2357  (1 
Curiously  enough,  that  paper  considers  a  f: 
(c)  with  varying  Fc££,  but  dot's  not  take  uj 
ary  case,  although  the  formalism  for  doing 
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As  justification  for  this  argument,  Figure  52  shows  an 
empty- resonator  calculation  of  I  >  |  vs  tilt  in  units  of  2  r 
for  F^jr  =  9.378  and  H  =  2.9.  The  validity  of  this  way  of  look 
ing  at  tilted  resonators  can  also  he  discerned  from  an  examina¬ 
tion  of  the  mode  patterns  themselves.  In  Figure  53  wo  compare 
an  exact  calculation  of  P(!2j-  =  9.3,  M  =  1.9,  2  c  Feff  =  1  with 
a  composite  formed  by  the  left  half  of  a  mode  pattern  for  T'{ ^ j- 
8.3,  M  =  1.9,  .r  =  0  juxtaposed  with  the  right  half  of  a  mode- 
pattern  for  Ft.fj-  =  10.3,  M  =  1.9,  c  =  0.  The  resemblance  is 
quite  closed  indeed. 
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Figure  53.  Comparison  of  Tilted  Resonator  Mode  with  Composite 
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SECTION  VI 


STABILITY  ANALYSIS  OF  A  SUPERSONIC  CHEMICAL  LASER 


A.  INTRODUCTION 

The  objective  of  this  study  is  to  develop  an  analytic  model 
of  the  stability  characteristics  of  a  supersonic  laser;  in  parti¬ 
cular  that  of  a  chemical  laser.  The  flowfield  and  optical  inhomo¬ 
geneity  of  the  chemical  laser  make  it  an  extremely  difficult 
problem  to  analyze,  and  it  becomes  necessary  to  simplify  the  laser 
characteristics  in  order  to  obtain  a  first-order  model  of  mode- 
medium  interaction  inside  the  cavity. 

The  chemical  laser  generically  consists  of  an  array  of  flow 
nozzles  from  which  alternate  layers  of  fuel  and  oxidizer  are 
ejected  (Figure  54).  This  loads  to  a  series  of  mixing  and  reac¬ 
tion  zones  across  the  laser  cavity.  The  subsonic  laser,  as 
described  previously,  can  couple  to  transverse  acoustic  waves  and 
result  in  instability.  The  supersonic  laser  will  generally  not 
have  such  a  mechanism,  as  the  lateral  transit  time,  D/cs,  is 
much  longer  than  the  axial  flow  time  l/'u,  such  that  acoustic  dis¬ 
turbances  will  be  swept  out  of  the  cavity  (see  the  glossary  of 
terms  at  the  end  of  this  section) .  For  the  purposes  of  this  work 
the  flowfield  is  treated  as  being  one  dimensional,  and  therefor 
any  coupling  of  disturbances  in  the  flowfield  to  the  lasing  wil 
be  primarily  axial  in  form.  This  ignores  the  possible  instability 
mechanisms  introduced  by  the  series  of  planar  mixing  layers  acting 
as  diffraction  gratings  or  prisms  which  could  then  interact  with 
the  laser  flux.  This  mechanism  is  considered  to  be  outside  the 
scope  of  the  present  analysis. 

B.  MODEL  ASSUMPTIONS 

As  stated  above,  our  model  treats  the  laser  as  If  it  has 
a  one-dimensional  flowfield.  This  is  obviously  not  true,  but  it 
leads  to  an  overall  simplification  of  the  problem  that  may  none¬ 
theless  be  reasonably  valid.  The  consequence  of  mixing  is  to 
reduce  the  overall  rate  of  chemical  reaction  in  the  laser.  Thus, 
a  reaction  of  the  form 


is  described  by 


A  +  B  -*  C 


dX 

-df  =  Wb 


where  X  denotes  mass  fraction,  and  p  is  gas  density. 


\  3-  1 
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Figure  54.  Schematic  of  Laser  Gcomotr 


Under  turbulent  mixinq  condition*',  the  local  values  of  the 
concentrations  [A)  and  [D]  may  vary  substantia  1 1 y  with  time,  and 
their  product  more  so.  Although  the  dynamics  of  the  mixing  pro¬ 
cess  are  very  complicated,  and  many  workers  have  investigated 
such  flowfields,  it  would  appear  feasible  to  use  a  simple  model 
to  estimate  the  reaction  rates  for  a  quasi  one-dimensional  treat¬ 
ment.  In  the  model,  the  concentrations  X;\,  XB  of  species  A  and 
B  (averaged  over  a  plane  normal  to  the  flow  direction)  are  given, 
by  their  mass-averaged  flow  rates  and  an  effective  reaction  rate 
as  follows: 


dX 
_ c 

dt 


VxaV 


where  X  is  a  mixing  parameter  which  is  typically  much  loss  than 
unity,  and  kfX  is  the  effective  reaction  rate.  The  rate  of 
disappearance  of  f uel /ox  id izer  species  can  be  related  to  a  lasing 
medium  length  in  the  flow  direction.  Thus  for  the  C7  laser,  with 
the  D2  concentration  significantly  greater  than  F2 ,  we  can  write 

X„  =  X  exp  (-k^pXX  t) 
r  i’  J-  L)  0 

O  2,0 

where  Xp  is  the  initial  fluorine  concentration.  Thus,  we  have 
an  e-folding  time  for  disappearance  of  fuel  given  by 


T 


kf pXXD 


2,o 


The  convective  distance  £  traversed  during  this  time  is 
thus  given  by 


£  =  UT 


where  u  is  the  convective  gas  velocity.  For  a  given  mixing  rate 
we  can  thus  relate  the  length  of  the  lasing  mccium  in  the  flow 
direction  to  the  kinetic  rate  and  the  flow  velocity,  such  that 


u _ 

kf  ipXD 


2,o 


In  the  following  we  denote  the  effective  reaction  rate  by 
kf  =  kf  X . 

A  chemically  ptimped  DF  laser  exhibits  multiline  operation 
such  that  laser  action  can  occur  on  an  many  as  three  vibrational 
transitions  and  a  total  of  15-20  vibra t  i  ona 1 - rota t i ona 1  lines.  (12) 


32.  Grom; ,  R.W.F.  and  Spencer,  D.J.,  "CW  Uydrogen-Hal  ide  Lasers.," 
Handbook  o_f  Chemic__a_l_T.ar._ers,  R.W.F.  Gross  and  J.F.  Dott,  ed .  , 
Wi  ley- In  terse  i  once  ,  New  York,  NY,  197G,  Chapter  •} . 
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In  order  to  accurately  model  such  a  system,  the  time  evolution 
of  all  lasing  states  and  the  flux  for  each  lasing  transition 
should  be  followed.  This  leads  to  a  system  of  approximately 
20-30  coupled  differential  equations  which  must  be  integrated 
numerically.  Such  a  large  set  of  equations  is  too  cumbersome  to 
be  of  use  for  the  present  application;  therefore,  simplifying 
assumptions  must  be  made  in  order  to  provide  a  manageable  kine¬ 
tics  scheme. 

The  first  assumption  that  we  make  is  that  three  vibrational 
levels  are  sufficient  to  model  the  behavior  of  the  lasing  species 
These  levels  are  the  upper  lasing  level,  DFU,  the  lower  lasing 
level,  DFi,  and  the  vibrational  ground  state,  DFg.  We  assume 
each  level  is  in  rotational  thermal  equilibrium. 

In  order  to  model  the  lasing  action  in  the  system,  it  is 
assumed  that  the  multiline  operation  of  the  DF  laser  can  be  de¬ 
scribed  by  a  single  transition  from  DFU  to  DF \ .  Implicit  in  this 
assumption  is  the  additional  assumption  that  no  lasing  occurs  to 
the  ground  state  of  DF.  This  allows  us  to  use  a  steady-state 
approximation  for  describing  the  unperturbed  lasing  medium,  since 
with  this  assumption,  the  ground  state  of  DF,  whose  population 
increases  with  time,  is  decoupled  from  the  remaining  kinetic 
equations.  Since  only  ~  20%  of  the  laser  enerqy  in  a  DF  laser 
arises  from  the  v  =  1  -*•  v  =  0  transition ,(  32  ,  3  3)  this  assumption 
is  not  unreasonable. 

The  reactions  necessary  to  describe  the  DF  laser  are  as 
follows : 


(i) Formation: 


F  +  D2  -*  DFu  +  D 


DF]L  +  D 

No  ground-state  DF  is  formed  by  this  reaction, 
(ii) Deactivation: 

kD 

DF  +  M  -*-U  DF,  +  M 
u  1 


D1 

DF,  +  M  -*  DF  +  M 

1  g 


33.  "NACL  Program  Final  Report."  (U),  J,TM-2‘J2,  TRW, 
Classified  :  CONFIDENT  IAI. . 


July  R,  1970. 
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M  is  the  combined  density  of  DF,  D2 ,  D  and  F;  deactivation 
by  the  carrier  gas  is  assumed  negligible. 


( i i i ) Lasing: 


DFu  +  hv 


g 

-*■  DF1  +  2  hv 


Three-body  recombination  is  assumed  neglibible  and  only 
single-quantum  deactivations  are  considered.  Estimates  for  the 
formation  and  deactivation-rate  constants  wore  obtained  from  the 
recommended  values  of  Cohen (34)  by  assuminu  that  the  D2/F  ratio 
was  2-3  and  20-30'i  of  the  F  atoms  were  reacted  in  the  lasing 
region.  Our  estimated  rate  constants  and  the  enthalpies  for  the 
various  reactions  are  listed  in  Table  1. 


One  more  quantity  is  necessary  to  model  the  behavior  of 
the  lasing  medium  and  that  is  the  gain  per  unit  length,  which 
can  be  written 


where  cs  is  the  stimulated  omission  cross  section  and  ANV/j  is 
the  population  inversion  for  the  v,j-l  -►  v- 1 ,  j  transition. 
Writing  this  equation  in  terms  of  the  species  discussed  above 
yields 


g  =  0  (  [DFU]  -  0  [DFJ_  ]  ) 

where  the  brackets  denote  species  concentration  and  0  is  the 
product  of  the  stimulated  emission  cross  section  and  the  frac¬ 
tion  of  the  total  DFU  population  which  is  in  the  lasing  rotation 
al  state.  For  a  P-branch  transition,  ••  -  exp  (-2J/Q),  where  Q 
is  the  rotational  partition  function.  Our  estimated  values  for 
0  and  0,  obtained  using  the  spectroscopic  data  contained  in  Ref. 
35,  also  are  contained  in  Table  1. 

Using  the  model  discussed  above,  the  time  evolution  of  the 
chemical  species  and  laser  flux  in  a  cw  DF  laser  is  described 
by  the  following  equations: 


34.  Cohen,  "A  Brief  Review  of  Rate  Coefficients  for  Reactions 
in  the  D2-F2  Chemical  System,"  TR-007 4 ( 45 30 ) - 9 ,  Aerospace 
Corp. ,  January  1974. 

35.  Emanuel,  G. ,  "Numerical  Model  inn  of  Chemical  Lasers," 

Handbook  of  Chemical _ l,a_se r_s ,  op.  cit..  Chapter  B. 
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TABLE  1.  CONSTANTS  FOR  SIMPLIFIED  DF  MODEL 


Constant 

kFu 

kFl 

kDu 


Best  Estimate 


0  .  n- 12  3  / 

8.3  x  10  cm  /sec 

9.4  x  10  ^cmVsec 

1.5  x  10  ^cmVscc 


Probable  Range  Heat  of 

5. 9xl0_12-l. 2xl0-11  -7.3 

6.6xl0_13-l. 3xl0-12  -23.4 

3.0xl0_13-4.5xl0”12  -8.0 


Dl 


8.6  x  10  ^3cm3/sec 


1.7xl0"13-2.6xl0_12 


-8.3 


o 


5.0  x  10_17cm2 


2.5xl0_18-2.5xl0-16 


0  0.5 


0.4  -  0.6 


Keaction 


kcal/mole 


kcal/mole 


kcal/mole 


kcal/mole 
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d  t  DF  ]  . 

- H_  =  k  l  F  ]  [  D  2  ]  -  kD  [DFuUM]  -  UDFu]  -  0  IDF*1)  ) 

d  t  u  u 


d  [°F  1  ntS 

k^lFjfD^  4-  k  IDFJCMI  -  k^lDF.JlM]  +  <  [DFJ -0  [DKj  ]  ) 


d  [DF  ] 

- 2-  =  k  [DF,  ]  [M] 

dt  1 


=  CO([DFu]  -  Q  [ DF x  ]  )  <J> 

where  c  is  the  speed  of  light,  hv  is  the  laser  photon  energy, 
and  $>  is  the  intracavity  flux. 

C.  ANALYTIC  MODEL 

The  laser  cavity  may  be  described  by  a  set  of  ono-dir<ions  i  on  a  ) 
flow  equations  plus  gain  and  flux  equation.  The  flowfiold  is 
assumed  to  be  inviscid.  Thus  the  continuity  equations  are: 


Continuity  3_0  +  3  (ou)  „  ~ 
3t  3x 


(212) 


Momentum 


3  u  ,  3  u 

P  vr  +  PUvm 


(213) 


Energy 


B£  -  [1W  KF1AIV1l''>(H2xr  +  ’W'AA,  ei-t) 


where 


+  K_  AH_  pX  X. 

ml 


3  ^  3 

3t  3x 


A2? 
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In  addition,  the  gases  arc  assumed  perfect  giving 

P  =  PRT  (215) 

The  reacting  species  may  be  represented  by  the  three  DF 
states  to  give  the  species  continuity  equations 


Upper  State 

DX 

u 

Dt 

v*?  V  Vv»  ’  ^  lx° 

-  OX,) 

(216) 

Lower  State 

DX, 

Dt 

1!F1DXFXH2+  T/Vu  +  R-;  lxu  - 

f‘Xl> 

(217) 

- 

K„  oxjC 

D,  ml 

Ground  State 

DX 

-  .2.  = 

Dt 

K  pX  X. 

D .  ml 

(21  8) 

The  lasing  may  be  described  by  the  flux  coupling  equation 
described  in  Section  II  and  the  gain  equation,  such  that 


Flux 


Gain 


3A 

3t 


CJ 

2L 


[ /u  u  '  dx  +  /u  T  g  '  dx  ] 


0  0 


y  =  (X  -  OX  ) 

DF 


(219) 

(220) 


by 


The  magnitudes  of  the  coefficients  in  Eg.  (219)  are  given 


a 


R 


2nd 


_ P _ 

PATM 


L 

A 


a  =  1.0 

These  equations  fully  describe  tlv-  f  1  uw  inside  t  he  cavity. 
For  the  purposes  of  a  stability  analysis  v:<-  are  ml  ere;  t . -d  in 
the  perturbation  of  the  flow  from  a  steady-st  ate  cond.  j  (ion.  Thus, 
re-wri  ting  each  varinbl  e_a  s  the  sum  of  a  steady  st  it.’  and  a  fluc¬ 
tuating  part  (  i  .  o.  ,  u  -  u  +  u'  )  wo  can  derive  p-  i  t  art  a  *  i  i -n  equ  t- 

tions  from  the  above,  non  loot  inq  hieti.  r  ordm  t  .  i  ns  in  flndint  in 
quantities.  For  t  hr'  purposes  of  this  analysis  it  becomes  :r  c<  : . . .  a 
to  consider  the'  mean-value  st  e.,dy- st  a  t  ,  prop,  a  t  -  .  <<;  tic-  f  b  >w 

variables  to  be  constant  through out  the  cavity.  There  is  cue 
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exception  to  this  case;  the  ground  state  mass  fraction  obviously 
increases  with  time  whereas  the  other  equations  have  both  source 
and  sink  terms.  The  set  of  perturbation  equations  becomes: 


Cont in  ui ty 


Bp'  —  3 o ’  x  -  3u' 

at  +  u  Tx  +  p  3^ 


(221) 


Momentum 


-  Bu’  ,  —  Bu'  Bp' 

p  at  +  pu  ax  =  "  ax 


(222) 


Energy 

Fx’+  §l’  =  "cL  [Kfuahfu+  kf1ahf11{xho-  l(V  V  xg)HxFo~(V  V  x 


(Kp  A I  Ip  +  Kp  AH  ]  { X.,  +  j  Xp  "  (X  +  X,+  X 
p  u  u  *1-1  Ho  1  ro  u  1  CJ 


!  J  ‘xu  4  xi  + 


+  ^  l\*\+  VVx>Xp  4  pxi 4  p'<v 


(223) 


where  XnQ  and  XpQ  are  the  initial  mass  fractions  of  oxidizer  and 
fuel  respectively. 


State 


P' 


pH  T'  + 
s 


P ' 


(224) 
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Species 


Upper 

DX' 
_ u 

Dt  ' 


\,<xH  -  2<V  V  5o,,,XF  ‘‘V  x1  +  V’0' 

n  ^  O 


-  KF  e(x  +  i  XF  -<x  +  V  V,,XV  XV  XV 

u  o  O 

DX’  , 


Dt 


Lower 


DX* 


DtT  -  Kf1  (XI1  -  2<XU+  V  Xg>  >  1XF  -<XU+  Xl+  Xg  1  >  p  ’ 

o  o 


-  K  5(X  4  J  X  -(X  +  X1+  x  )(X'U+  X'1+  X'  > 

loo 


+  KDu<V  Xm)  P '  +  KDuD  Xm  X '  u 


+  &  i<v  i<**u-5*v 


'd/Vi"'-  KD1<’Xb.  X'l 


Ground 


DX'  „  9X 

——2  =  KD,  [X-  X,o’  +  X  p  X'  ]-  u '  FT  -2 
Dt  lml  m  1  dx 


Gain 


(  (Xu-  OXjJp*  +  p X ' u  -  fOX'j  < 


pX^O  In*’ 


T  '  ] 


DF 


I/O 
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(22! 


(22 


(22 


(2? 


r  ,s  i  i 


1 


Flux 


9 1 


2L 


[  r  aRp  '  dx 
o 


/  ctjg'dx] 


(229) 


As  wo  are  considering  supersonic  flow,  we  need  to  only  con¬ 
sider  the  flowfield  downstream  of  the  nozzle  bank;  thus  a  Laplace 
Transform  treatment  appears  tractable.  We  apply  a  double  trans¬ 
form  in  both  space  and  time;  i.e., 

2  ,  ,  .  -oo.oo  -st  -kx,,,  .  .  .. 

f(s,k)  ~  f  j  e  e  f(t,x)dxdt 


This  produces  a  set  of  linear  equations  from  the  above 
partial  differential  ones.  Thus  Eqs.  (221)- (229)  give 


Continuity 


Momentum 


(s+uk)  p  +  p  ku  =  0 


p  (s  +  uk) u  +  kp  -  0 


(230) 

(231) 


Energy 


1  ,  r. 


-UKp  AHd  +  Kr  AH,,  HX„  -  j  (Xu+  X1+  Xg)}(xF  -  (Xu+  Xg)} 


F  F  F,  F.  H 
u  u  11  o 


+  [K„  (  L  +  K  AH  ]  x  +  X,]X  }—  p 
D  D  D.  D,  U  1  m  c 
u  u  1  1  P 


+  (Ikf  a“f  +kf/"f ,Hxh  +  2xf  "  'V  V  V>-Ikd  ‘"h'W*.1 

uull  o  o  uull 


x  —  {X  +  X.  } 
c  u  1 
p 


(232) 


+  (s  +  uk)T  - k  u/c^p  p 


*  5T"V  A"f  +  *y*»r,  >  *  IXF  -«V  V  V’  *« 

puullo  O 


=  0 
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State 


—  p+pR  T  -  p  =  0 

y  g 


(233) 


Species 


Upper 

State 


1,5 


-  i'V  xl  +  y)lXF  '<V  V  V1* 

u  o  o 


1  X. 


+  [  (s+uk)  +  Kf  p(XF  +  2  F0-<V  Xl+  xg)nx 


F  "  “*F 
o  o 


(234) 


+  (  (s+uk )  +  Kp  p(XF  +  |xF  -  (Xu+  Xx+  Xg) )]X1 


u  o 


o 


+  kf  o  (x  f  x  -(f:u+  x1+  x  )>x 
u  o  o 


g  g 


Lower 

State 


3<V  V  WlV  "1+  V’*  \V»-  V*! 


+  (K  P(x  +  2xf"(Xu+  Xl  +  Xq))  “  Kn  r'Xm  hv  ^  Xu 
F1  Fo  o  y  DU 


H^V  °X1>  + 


(23S: 


+  [KF  p(x  +  |  XF  -(Xu+  Xx+  X  ))  +  (s+uk)  +  h-~  +  kd  ;xmi 

1  o  .o  1 


X  1  X. 


♦'£  T4' ; 


♦  x  MX,,  *  2  XF  -'V  V  V’  X9 
1  o  o 


(i  round 
St  ate 


-(s+2uk)KDiXniX1C  -  X.UK  <  uk  ( r  )  uk )  Xy  »  0.  (2  3C 
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r 

« 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 


Ga  in 


g  -  M—  (<XU  -  0X: )  c>  +  P  Xu  -  P  v>:x  + 


(  X  L  In' 


T]  “  0  (23  7 ) 


Flux 


sC  -  ^  {  -»'R  o  +  a, g }/k  =  0 


(238) 


F.1  imi  nat  ing  u  from  Eqs.  (230)  and  (231)  and  c;  from  Kqs. 
(237)  and  (238)  wo  obtain 


Continuity-Momentum 


(s  +  uk)^  p  -  k^p  =  0 


(239) 


Gain-Flux 


§L  (V  °i5(V  3*i)lp  +  T£  riip  xu  *  ks; 


-  cf  Z  ,  ~  c-’  a  t  r’  X 1 c'  Ire-  n 

o  2L  ctjpt  X1  +  o  _I. — _ -  J  0 


(240) 


whore  a  =  a/M  . 

The  seven  equations  in  variables,  p,  T,  p,  <,  X  ,  X^,  Xfi 
can  be  written  in  matrix  form 


[A]  (x]  =  0 


(241) 
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Thus  wc  have 


Trans  foiTieO  Kg  nations 


a 

11 

0 

0 

0 

a 

15 

a21 

a22 

-ks 

Ct24 

a,r 

(s+uk ) 

2  0 

0 

0 

0 

a41 

a42 

0 

a  .  . 

4  4 

(s+uk) 

a5l 

a52 

a5  3 

015  ^  +  ( s+uk) 

n5  5 

a61 

a,. _+ (s+uk) 

0  L. 

0 

ag4+ (s+uk) 

0 

a71(s+2uk)  ° 

0 

a74^ 

0 

0  0 


-k  0 


a  *  /-  k  a  — * 

46  47 


0  a6  7 


0  uk (sxuk) 


(242) 


where 


11  =  a15  ~  pRg '  a16  "X 


_  £i 

21  2L  1  ~R  "1“  '"u 


u,  =  —  [a„+  aTo(X„-  exx)1 ,  a 


=  a  r 

22  2L  uI! 


24 


—  c<J>  -  pxl01n0 

2L  ai6pa  '  a25  =  “i  2L  “  ~~  ~ 


T 
1  ,= 


41 


=  -  ^-UK  Ml  +  K  MlF  HXh  -  2'(XU+  Xx+  Xg)}(X  -(Xu+  Xg)} 

v'p  U  U  1  1  o  ° 


+  (KD  A!!D  +  KD,AI,D,HV  V  VXm  ] 

u  u  11 
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42 


a44  =  c~  ^ 1K^  A” 


F  ""F  ”1  -i 

u  u  i 


+  K  AH  Hx  +  '2  X  ~(XU+  XL+  xg)} 

I  i  i  r\ 


1  “o 


•V!V\“Srs'' 


46 


u 

pc 


=  JL{y  +  \  X  -(X  +  X,+  x  )  }  (Kj,  .Hj,  \\]  ] 

47  C  1  H0  2  Eq  U  1  q  u  u  1  1 


P  47  Cp  0 


-IK  (XH  "  \  (V  Xl  +  Sg))(XF  _au+  Xl+  Xg)  +  KD  XmXu 


51 


1  o 


a  2  =K  P(X  +  |  Xp  -(Xu+  Xx+  Xg)-K  P  Xm  hv 
J  1  o  o  u 


0  (x  -  ex,) 


L.  .  —  r —  ->  j 

53  hv  u  1 


V,  -*F  i«,  +  1  *F  -'V  V  V>  *  +  V  X” 

loo  x 


o^j  QlnO 
a55  hv 

a57  »Kf.  5(Xh  *  ixr  -(V  v  y> 

J  X1  o  o 


«tl  -  -xP  P  (X„  -  I  'V  V  S,"«r0-(V  V  V1 

u  u  o  ° 


•S2-“64-“67-Kru'P<V’V(5U+  V  V’ 


a71  =  “VXmX.  '  a'7  4  '  hD1Xmp 
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The  stability  of  the  system  can  be  examined  in  terms  of  the 
characteristic  equation  derived  from  Eq.  (241);  i.e., 


det[A}  -  0 


(24  3) 


An  analysis  stability  of  the  system  is  obtained  from  the 
Routh  criterion  for  Eq.  (  243).  (36) 

The  characteristic  equation  may  be  obtained  from  Eq.  (242) 
After  tedious  algebra  this  reduces  to 


a^s 

;7  +  a6 

6  .  5 

s  -r  a gS 

+  a4s4  + 

a3s3 

+ 

2 

a2s 

+  alsl 

+  no 

where 

u? 

A6 

a6  = 

6  qAg  + 

A5 

a5  = 

15n2Ar 

+  5r|Ag  + 

A4  +  b5 

a4  = 

20»i3Ag 

+  10q2A5 

+  4tiA4  + 

A3  + 

5nB5  + 

B4 

a3  = 

15r’4\ 

+  10i|3Ag 

+  6qA4  + 

3r,A3 

+ 

A2  + 

10m2B5 

+  4nB4 

a2  = 

6n5A6 

+  5q4Ag  + 

4n3A4  + 

3.,2A3 

+ 

2nA2 

+  A^  + 

i0n3B5 

+  6n  +  3nB3  +  b2 


6  5  4  3  2 

-  h  Ag  +  n  A5  +  n  a.  +  n  a3  +  n  'A  +  qA,  +  a 


0 


+  5n4B5  +  4n3B4  +  3n2B3  +  2qB2  +  13  3 


5  4  3  2 

-  n  Dj  +  n  +  n  +  n  1^2  +hb^  +  Bq 


36.  DiStofano,  J.J.,  Stubbcrud,  A.R.,  Williams,  I.J.,  "Feedback 
and  Control  Systems,"  McGraw-Hill  (1967). 
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* 

a 

t 

f 

a 

a 

a 

a 

■  ■ 

a 

a 

a 

a 

i 


and 


a6  =  -a16kn 


A5  =  a15a46k  n  "a16a62kn  +a16 (a52  "  a54)kn 


B5  a16a53  (a22  "  a24,n 


A4  allk  n  a15a46(a52  '  a54)k2n  +ui6a62  (a52  a54 

+a16a57a74kn 


B4  a53(a22  a24) n ^“a15a^6k  +  a621 


A3  allk  n^a52  a54^~a62^+0t15a46a57a74k  n 

+Ct15a47a71k  -  a15a41nk  +a16a57a62a74kr| 


a16a47a55a74kn  a16a52a62a74kn 


B3  alla53(a22  “a24)k  n  +  a15a46a53a62  (a22  "a24)kn 


a16a25a47a53a74n  a16a2 2n53aG2n7 4 n 


)kr, 
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alla62  (a52  "  a54)k  n  +  a57cx74k  r' 


+a15a46a57a62a74k2n  +  a15a4 7a62a71k3 


~a15a47a71 (a52  “  a54)k  +  a15%7a71k  n  _a15a42a62 


15  47  7 1 


+a15(a41a52  "a41a54)k3n  +  a^a4fia52a62a74k2’' 


1546526274 


-a15a42a62a71k2  +  a15a42a61k3n 


+a16a55a62a74la42  “  a47)kn 


alla53a62(a22  ~  a24)k  n  a15a47a53a71  (a22  ~Cl24' 


+a15a41a53  (a22  "  a24)k  n  +  a15a22a46a53a74kn 


+a16a25a53a62a74 (a42  a47)n 


alla47CX55a74k3n  ”  alla52a62a74k  n 


+ct15a47a62a71k3n  r  a15a47a62a71  (a52  "  a54)k 


-a15a47a71(a52  "  °54)k  n  +  °l5a7 4 (a4 1°5 7n5 1*4 7 ]  k 


+a15a62a41(a52  “  «54)k3n  "  a15a42ftfa2Cl71k3n 


+a15a42°62a71  (n52  ~  a54)k3  ~  a15a42rt61  (a52rt54  }  k 


1  11  25  47  53  74  1  11  22  53  62  74 

“ai50t47a53Ct71  (a22  "  a24)k  n  "  (t15a47Cl53C,G2a71  (ct22  ~  a24)k 

“a15a21a4'a53a74k  n  +  a15a41a53a62  (a22  “  az4)K  n 
+a15Ct42a53a61  (a22  "  a24)k  n  +  a15a42a53a62a71 (a22  ~  a24)k" 


0  =  "alla^7a55aG2a74k3n  ~  a15a4 7a62a71  (a52  "  a54)k3n 


+Ct15a62a74  (n41a57  “  a51a47)k  n 


+alla42a55a62a74k  n  a15a61a74  (a42a57  K47a52)k  n 


“a15a62a74 {a41a52  “  a51a42)k  n 


+a15a42aG2a71 (a52  ~  a54)k  n 


-alla25a470f53aG2a74k  n  “  ai5a21a47a53aG2a74k  n 
-C115a47a53aG2a71  (a22  '  a24)k  n  +  alla2  5a4  2a5  3aG  2n74k  n 

+Ct15Ct53a74  [a22a47a61  +  a62{ft21U4 2  "a220l4 1^  1  k  n 


+Cl15a42a53a62a71  (a22  “  a24)k  1,1 


The  va]ua  of  the  spatial  wavonumbor ,  k,  is  taken  to  be 
given  by  the  cavity  dimension  in  the  flow  direction,  •.  Thus, 

k  1/1.  ...  . 

Thus,  v.o  can  now  numerically  specify  all  the  reef  f  ic  i  ent.s  of  tin- 
character  i  s t ic  equation. 
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D.  STABILITY  ANALYSIS  OF  DF  LA SLR  MODEL 


The  values  of  the  laser  parameters  can  be  obtained  from 
the  "gain  equals  loss"  condition  and  Eqs.  (21G)  and  (220).  The 
"gain  equals  loss"  condition  specifies  the  gain  in  terms  of  cav¬ 
ity  loss  and  cavity  length  such  that 


e2gL 


r2 


where  rj  and  i'2  are  the  mirror  reflectivities.  Thus  we  obtain 
the  population  density  of  upper  and  lower  states  from 

(Xu  -  0X1)  =  g/cp 


The  flux  can  be  obtained  from  Eq.  (216)  where  we  assume 


DX 

u 

Dt 


=  0 


0X1) 


Thus  4>  =  K  P2X  X  ahv/(og) 

Fu  Fo  D2,0 

Using  these  criteria,  the  parameters  for  a  wide  range  of 
flow  conditions  were  investigated.  The  laser  conditions  are 
listed  in  Table  2.  Over  the  entire  range  of  mixing  rates,  flow 
velocities,  cavity  lengths  and  pressures,  use  of  the  Routh  cri¬ 
terion  specifies  that  the  laser  is  stable.  This  can  be  simply 
described  by  reference  to  the  coefficients  of  Eq .  (24*1).  The 

only  matrix  coefficient  containing  contribution  from  (the 
density  fluctuation  induced  flux  variation)  is  '*21  which  appears 
in  coefficients  and  Bq;  these  affect  only  coefficients  a]  and 
aQ.  In  general,  Bi  and  Bg  are  significantly  smaller  than  the 
other  terms  in  ai  and  ag,  resulting  in  both  ai  and  ag  being  posi¬ 
tive  (if  either  of  these  coefficients  were  negative  an  unstable 
root  would  occur) . 

The  response  of  the  laser  to  upstream- induced  noise  dis¬ 
turbances  may  be  described  in  terms  of  transfer  functions  derived 

from 


[A] 

lx]  - 

IF] 

where 

[F]  is  the 

for  c 

i  ng 

f  uricti  on 

vector . 

The  transfer  function 

between  parameter 

xn 

and 

Fn 

is  given  by 

y-n 

dot. 

(A,,] 

F- 

dot 

Fa  P 

where 

( A  p ]  is  the 

ma  1 

r  i  x 

[A] 

wi  th 

the  nfh 

column  replaced  by  (F] 
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TABLE  2. 


LASER  CONDITIONS  INVESTIGATED 


Cavity  length  (L)  0.1  -  3.0  m 
Mixing  length  £  1.0  -  10.0  cnis 
Flow  velocity  10^  -  3  x  10^  m/sec 
Cavity  pressure  1-10  torr 
Cavity  density  5x10^-5x10^  kg/m^ 
Laser  mix  10  He/5  D2/O.8  F2 
M  =  5.45 
Y  =  1.59 

Mirror  Reflectivities 

r  ^  =  0.98  r2  -  0.88 
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The  response  of  such  a  .system  is;  thus  characterized  hy 
1  /dot  (A)  .  The  normalized  inuq.n  i  t.udcs  of  this  f  uj.ct  i.  on  .in1  i  1  1  u:; 
t rated  in  Figure  D5.  The  response  is;  essentially  uniform  up  to 
frequencies  of  105  rod /sec ,  and  thereat  ter  rapidly  dec  rease:;  in. 
magnitude  with  a  slope  a svinp toting  to  -■.*'.  Tiie  response  reduce 
in  magnitude  as  cavity  pressure  is  increased. 

Obviously,  this  is  only  part  of  the*  overall  transfer  f unc 
tion,  and  further  analysis  is  necessary  if  the  detailed  system 
response  is  to  be  determined. 

E.  DETAILED  SYSTEM  RESPONSE 

The  model  developed  here  is  a  simplification  of  the  real 
laser  behavior.  A  more  detailed  and  correspond inc] y  difficult 
approach  is  required  if  the  full  threo-dimonsjonu 1  nature  of  th 
mode  medium  interaction  is  to  be  investigated.  From  t he  prev i o 
analysis  a  simple  acoustic  or  gain  f  luctuation  inst:  ibi  1  ity  does 
not  appear  to  occur  in  the  chemical,  laser.  However,  the  real 
flowfield  may  result  in  other  significant  mechanisms  occuring. 

The  chemical  laser  consists  of  a  series  of  quasi-; lanar 
flame  sheets  (Figure  54).  The  jets  from  the  nozzle  bunKk  are 
normally  in  the  transition  range  of  Laminar  to  turbulent  flow, 
’..'here  laminar  instabilities  of  the  f  ]  amoshoof  are  apparent. 
These  instabilities  arc  akin  to  the  Ray leigh-Tay ! or  instability 
in  mixing  layers  (see  for  example,  Kimura(37)  acd  'l’oonq .  ( 3  8 ) 
Fluctuations  in  the  local  gain  and  flow  proper  t  i  or,  will  result 
in  modulation  of  laser  output  flux.  Whether  such  a  f 1 ow field 
instablity  can  lead  to  an  output  instability  is  not  clear.  The 
coupling  between  the  array  of  mixing  layers  is  of  importance; 
if  the  flowfield  is  driven  by  an  upstream  acoustic  disturbance 
which  is  identical  for  all  the  shear  layers,  this  may  result  in 
significant  optical  coupling.  However ,  a  random  phasing  of 
oscillation  of  the  array  of  flame  sheets  may  result  in  a  quite 
different  behavior. 

F.  SUMMARY 

A  simplified  model  of  a  supersonic;  chemical  laser  he  be.* 
derived  and  a  stability  analysis,  performed.  This  indie. ,t  *  t  ha 
the  laser  is  stable  -across  a  wide  range  of  operating  eond  i  t  ions 


37.  Kimura,  I.,  T<*nt  h  Symposium  ( Int  .  )  or.  Combustion,  p.  1  — 

1  300  ,  The  Combust  ion  Institute  (1°65). 

38.  Toong,  T.Y.,  Sulont,  R.F.,  Stanford,  J.M.,  AruL  i  ..on  ,  c.y., 
Tenth  Symposium  (Int.)  cm  Combust  ion,  p.  1301-1  Hi,  Th** 
Combustion  lns.fi  tut.  (lOf.S). 


AVCO  ! 


However,  the  simplifying  assumption:;  r<*niov->  the  po;u;  i  h  j  1  itv  of 
a  coupling  occurring  an  a  consequence  of  f  1  nine- sheet  instabi  1  i  ty 
and  as  such  would  not  appear  in  the  solutions  obtained  in  this 
work . 


The  above  analysis  provides  a  useful  starting  noint  in 
obtaining  a  description  of  possible  mode-mod  i  uni  interactions  in 
supersonic  mixing  lasers.  However,  it  in.  far  from  complete  and 
it  is  recommended  that  further  exploration  of  the  possible  be¬ 
havior  of  flame-sheet  instabilities  upon  optical  performance  is 
necessary . 


G.  GLOSSARY  OF  TERMS 
a  coefficients  in  Eq.  (244) 

A  coefficients  in  Eq.  (244) 

B  coefficients  in  Eq .  (244) 

c  speed  of  light 

cg  sound  speed 

Cp  specific  heat 

g  gain 

n  enthalpy 

k  space  transform  variable 

k  kinetic  rate  coefficient 

K  effective  rate  coefficient  (=  Ak) 

1  mixing  length  scale 

L  cavity  length 

Hnr  molecular  mass 

r  mirror  reflectivities 

p  gas  pressure 


u  gas  velocity 

t  time 
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Ill 


s 


time  transform  variable 


T  temperature 

x  axial  coordinate 


0tl'flR  flux  coupling  coefficient  in  Eq .  (219) 

a  coefficient  in  matrix  Eq .  (242) 

3  Gladstone  Dale  Constant 

Y  ratio  of  specific  heats 

0  degeneracy  coefficient 

X  mixing  parameter 

X  laser  wavelength 

p  gas  density 

PATM  density  at  a  pressure  of  1  atm 
o  lasing  cross  section 

<t>  intracavity  flux 

t  reaction  time  constant 

X  mass  fraction 

AH  enthalpy  of  reaction 

o  o/Mdf 


Subscripts 
u  upper  state 

1  lower  state 

g  ground  state 

F  formation  rate 

D  deactivation  rati' 
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SECTION  VII 

SUGGESTIONS  FOR  FUTURE  RESEARCH 


The  dynamic  models  considered  in  Sections  II,  IV  and  VI  have 
made  simplifying  assumptions  on  the  nature  of  the  Inner  cavity. 

In  particular,  uniform  values  of  gain,  densit  y  and  f  1  u>:  have-  L-e*.  :s 
assumed  in  the  cavity.  More  quantitative  results,  may  be  obtained 
by  perturbing  the  actual  steady-state  modes. 

The  theory  contained  in  Sections  Il-IV,  and  VI  of  tins  re¬ 
port  deals  not  only  with  possible  system  instabilities. ,  but  also 
in  the  system  transfer  functions.  Accordingly  statistics  of  the 
system  output  can  be  predicted  if  external  perturbations  can  be 
modeled.  In  particular,  it  would  be  useful  to  model  the  turbu¬ 
lent  nature  of  the  nonuniform  flow  field  of  the  flowing  gas. 

Once  the  correlation  (or  distribution)  function  of  the  density 
(i.e.,  refractive  index)  function  is  known,  the  stn  t  i  sj  t  ion  of 
tho  noise  driver  are  known.  The  theory  o:  this  report  then 
allows ,  via  the  transfer  function,  the  calculation  of  the  sta¬ 
tistics  of  the  laser  flux  output.  In  addition,  tho  net  hods  of 
the  theory  of  light  propagation  through  turbulence  can  be  used 
to  find  the  static,  i  tics  of  wavo-f  rent  di  start  ion  of  the  out  rut 
wave.  This  problem  is  especially  acute  for  the  chemical  laser, 
where  the  imperfect,  mixing  of  the  reactant  gases  can  bo  modeled 
to  predict  the  statistics  of  the  system  inhomoyeneities . 


calculating  the  sell-consistent  mad  os  c:  a  saturate  a,  v.r.at:.;.  . 
resonator.  This  is  in  fact  the  only  reasonably  : ast  technique 
that  we  know  of  for  making  such  a  calculation.  Althouyn  we  have 
only  appl  ied  it  to  a  str i p- resonator  cs  come  try,  it  appears  to  be 
genera  1 i cable  to  a  circular  geometry  as  well.  Calculations  of 
this  type  can  bo  used,  for  example,  to  determine  the  coup! inn 
between  modes  and  the  modification  of  tho  out  put  coup1,  i  ng  of  in¬ 
dividual  modes  induced  by  an  externa]  perturb.! t  ion.  v.V  predi  ct 
that  this  technique  will  find  wide-ranging  appl ieat ions  to  laser 
resonator  problems  in  the  future. 
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APPENDIX  A 


THE  PERTURBED  RESONATOR 


In  this  appendix  we  will  treat  the  perturbation  expansion 
of  Section  II  somewhat  more  riqorously,  therobv  cla rifyino  the 
physical  meaning  of  the  acoustically  induced  perturbation  o£  t ho 
electromagnetic  wave  in  the  resonator. 

Writing  the  part  of  the  electromagnetic  wave  which  propagate.1 
to  the  right  as 


E(x,  t) 


u  ( x  ,  t )  e 


i(k0z"'ot> 


+  c.c 


(A- 1  ) 


we  find  that  u  satisfies 


u(x,t  +  i  )  =  f  K(x,x',t)  u  ( x 1  it-)  ax' 


(A- 2) 


where  K  =  K  +  K x ;  Kq  is  the  unperturbed  propagator  (i.c.  ,  t he 
propagator  in  the  absence  of  the  acoustic  wave!  ;  is  the  per¬ 
turbation  due  to  fluctuations  in  the  density  (acoustic  wave)  and 
gain  of  the  medium  (see  Section  II);  x  and  >: '  represent  (vector) 
positions  in  an  x-y  plane  between  the  mirrors;  and  :Q  =  2L/c  is 
^he  round-trip  time.  In  the  following  we  write  x  to  denote 
x  =  (x,y) . 


The  time  dependence  in  is  due  to  an 
whose  frequency  is  very  low  compared  to  l/~0 
means  that 


acoustic  oscillation 
-  5  x  107  sec-1.  Thi 


3K. 


9 1 


o 


(A- 3) 


and  we  therefore  may  construct  an  equation  for  3u/~t  (x,t)  as 
f ol lows : 


?u 

at 


~~  t  u ( X ,  t  +  I o ) 
;  o 


n  (x  ,  t) ]  = 


+  /  V  j  ( x  ,  x  '  ,  t ) 


u (x '  , t)  dx  '  ) 

a** 


i  l/Vo(x,x’)  u(x'.t)  dx' 

(A -4) 
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where 


V 

o 

s  “  (K0(X'X,)  " 

6  (x-x' ) 1 

(A-  5 ) 

To  ° 

and 

V1 s  f  Vx'x' 

o 

1  ,t) 

(A-6) 

Equation  (A-4) 

has  been  written 

in  a  form  which  makes 

i  t 

identical  with  the  Schrodinger  equation.  Thus,  we  may  solve  this 
equation  by  the  well  known  perturbation-theory  techniques  used  in 
quantum  mechanics. 

The  integral  operator  VD  represents  an  unperturbed,  time- 
independent  "Hamiltonian"  of  the  system,  while  is  tne  perturba¬ 
tion.  To  calculate  the  behavior  of  u  (>: ,  t)  (analogous  to  the  wave- 
function  of  quantum  mechanics)  w-e  may  expand  u(x,t)  '  any  complete 
set  {'in(x)}  of  space-dependent  functions: 


u  (x,t) 


an(t> 


V 

'  n 


(x) 


(A-  7 ) 


As  the  set  {yn}  we  choose  the  eigenfunctions  of  the  unperturbed 
Hamiltonian  VQ,  so  that 


AH'  (x)  =  /  V  (x ,  x  ’  )  T  (x')  dx'  (A- 8) 

nn  o  '  n 

i.e.,  H,n  is  the  eigenf unction  of  eigenvalue  \  .  Equation  (A-8) 
may  be  shown  to  be  equivalent  to 


t  (X)  =  (1  +  il  r  )  /  K  ( x , x 1 )  7  (x'j  dx’ 
n  n  o  o  n 


(A- 9) 


or,  for  |AnTo]<<  1 


-iA  r 

e  n  °y  (x) 
n 


/  K  ( x  ,  x  1  )  j'  (x’  )  dx' 
o  n 


(A- 10) 


m 

/S'/ 
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which  is  the  usual  integral  form  of  the  eigenvalue  equation  of 
a  laser  resonator.  It  will  be  shown  below  that  the  natural  optical 
frequency  of  mode  n  is 


u  =  uj  +  Re  A 
no  n 

i.e.,  Re  An  represents  the  frequency  offset  of  mode  n  from  that 
of  the  pure  longitudinal  mode  of  frequency  . 0 .  On  the  other  hand, 
-2  ImAn  represents  the  fractional  energy  loss  per  second  of  mode  n. 

Substituting  Eqs.  ( A—  7 )  and  (A-8)  into  Eq .  ( A—  4 )  gives 

‘Efe 

m 


l  A  a 

m  m 


T  (x) 
m 


E 


m 


a  (t) 
m 


/  V.  ( x  ,  x  1 


,t)im(x 


)  dx1 


We  now  assume  that  the  modes  { M'  n }  form  a  complete  ortho¬ 
normal  set.  This  assumption  could  be  rigorously  justified  if 
the  unperturbed  Hamiltonian  VQ  were  Hcrmitian;  the  eigenvalue 
An  would  then  be  real,  corresponding  to  an  unperturbed  cavity 
in  which  all  modes  have  the  same  value  of  net  loss  (loss  minus 
yuin).  This  will  clearly  not  be  the  situation  in  an  actual  un¬ 
stable  resonator.  However,  the  only  modes  which  arc  interesting 
are  the  initially  oscillating  mode  and  those  modes,  who.se  net  losses 
arc  nearly  equal  to  that  of  this  initial  mode;  other  modes  (with 
significantly  higher  net  losses)  will  bo  strongly  damped  and  will 
therefore  never  be  excited  by  the  acoustic  coupling.  Thus,  .it 
should  be  a  good  approximation  to  assume  that  for  the  modes  of 
interest 


1mA  -  0 
n 


(A-  1  2 ) 


and  therefore 


/  T  * 

n 


V  dx 
m 


=  6 

nm 


(A-  1  3  ) 


Equation  < A— 1 2 )  implies  that  multimode  oscillation  is 
model  laser.  In  fact,  such  multimode  behavior 


in  our 


observed  in  ew  Hill,  husc-rs 
by  T  *(x)  and  integrating 

Eq.  (A- 1 3 ) 


.  (if 

over 


Multiplying  both  sides 
all  x  then  yields,  with 


i  'OSS  i  b  1  e 
has  be>ui 
e!  Eg.  (A- 1  1  ) 
the  a  1 d  of 


1.  R.  Patrick,  AKRb,  private  communication. 


ir? 
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(A- 1 4  ) 


da 
_ n 

dt 

where 

Wnm(t)  =  ff  4’n*(x)  M'm  <x  ‘  )  dxdx’  (A-15) 


i  A  a 
n  n 


-iZw 

* —  nm 

m 


(t) 


a  ( t ) 
m 


It  will  be  noted  that  no  use  has  been  made  so  far  of  perturbation 
theory;  within  the  limits  of  our  assumption  that  the  modes  of  in¬ 
terest  are  orthonormal  Eq.  (A— 14)  is  exact . 


The  meaning  of  Eq .  (A-14)  may  be  elucidated  by  first  ex¬ 
amining  its  solution  in  the  absence  of  the  perturbation  (i.e., 
with  Wnm  =0).  We  then  find  that  each  amplitude  an  satisfies 


da 

— rp  =  -  iA  a 
dt  n  n 


(A- 16 ) 


with  the  solution 


-i A  t 

a  =  a  (0 )  e 
n  n 


(A- 1  7 ) 


so  that  the  electromagnetic  wave  propagating  to  the  right  in  the 
cavity  is,  from  Eqs.  ( A  —  1 )  and  (A-7) 


E(x,t) 


(0) 


e 


i (k  z-w  t ) 
o  n 


V*) 


+  c .  c 


(A- 1 8 ) 


where  w  ,  the  natural  frequency  of  mode  n,  is  given  by 


u>  =ca  +  I<e  ( A  )  =  o)  +  A  (A- 19) 

no  non 

This  i s  clearly  what  one  would  expect:  each  mode  which  is 
present  initially  (at  t  =  0)  continuer,  to  oscillate  freely  at  its 
natrual  frequency.  Becau.se  we  have  not  allowed  for  any  dumping, 
none  of  these  modes  decays  with  time;  since  we  have  not  allowed 
for  any  coupling  between  them,  they  do  not  affect  one  another, 
cither. 
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Consider  tho  direct,  effect  of  the  perturbation  on  an  init¬ 
ially  oscillating  mode.  We  denote  f he  initial  mode  by  the  sub¬ 
script  0,  and  take  \Q  —  0.  (This  is  equivalent  to  defining  ,() 
to  be  the  frequency  of  the  initial  mode.)  I'i  rst -order  :>.-rturb  :t 
theory  allows  us  to  ignore  the  amplitudes  of  all  the  other  modes 
to  a  first  approximation.  We  thus  have,  from  Eq .  (A— 14) 


da 

o 

dt 


i  W  (  L  ) 

oo 


a 

o 


(A- 20 ) 


The  real  part  of  N0q  (defined  in  lb:.  (A- 15)  )  correspond:; 
to  a  shift  in  the  frequency  of  the  initial  mod',  while  the  im¬ 
aginary  part  corresponds  to  a  change  in  tho  :iet  cavity  loss  (or 
gain)  per  unit  time  of  this  mode.  Sued;  a  loss  chan  .-••’could  core 
about  as  a  standing  transverse  a count ic  wave  osci list,  s  in  ampli¬ 
tude;  we  would  tiler,  expect  the  cavity  loss  to  be  modulated  at’ 
twice  the  frequency  of  the  acoustic  wave.  Loss  modulation  of  t his 
type  is,  in  fact,  used  to  mode  lock  laser s.  It  could  si  princn 
lead  to  an  instability  if  the  initial  mode  is  not  uniform  in  an*, 
sity  within  the  cavity,  since  the  heating  associated  with  t  1  ■  e.-.vity 
loss  modulation  is  then  nonuniform  an  1  can  create  an  acoustic  dis¬ 
turbance.  Indeed,  our  calculations  of  the  mod. -s  of  an  unstahm 
resonator  (Section  V)  show  explicitly  that  a  nenun  i  form  node  i 
sity  distribution  is  to  be  expected,  so  this  is  a  possible  candi¬ 
date  for  an  acoustic  instability. 


Let  us  now  consider  the  effect  of  the  acoustic  :  ert.urb  it  ion 
on  the  other  modes  of  the  cavity.  Again,  wo  may  assume  to  fitst 
order  that  the  amplitude  aQ  of  the  initial  node  is  much  luma 
tlian  that  of  any  other  mode.  We  may  thus  ignore  all  amplitudes 
am  other  than  a0  on  tho  right  of  Eq.  (A—  1 4 )  ,  and  set  e;u  >.l  to 
unity.  This  approximation  is,  of  course,  only  self  s  :ix :  » 

as  long  as  the  equations  predict  that  ' a  ;  ■<  1  fer  all  n  v  n, 
i.e.  ,  in  tho  initial  stages  of  any  buildup'  of  other  cavity  rod.  .  . 
The  aR ' s  then  satisfy 


da 


I 


n 

dt 


+ 


i.\  a 

n 


n 


i  W 


no 


(t)  aQ(t) 


(A- 21 ) 


It  is  useful  at  this,  point  to  consider  two  ty:  •  • 
which  we  label  K  and  M  for  def  ini  teness.  Ms  N  is  e 
frequency  with  the  initially  present  mode  (  0 )  ,  w 
is  offset  in  frequency  by  an  amount  M  f  0.  The  mod., 
satisfy 


daN 

dt 


-  1  vWl) 


o ;  i’i.  ■  i .  •  s , 
at  •  ■  , 
:).•  :m-i. 
umpl  i  tud.  o 


(A- 22) 
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Associated  with  ay  will  be  the  perturbation  in  t he  field, 
n.j  (t ) 't'N  (x )  ,  which  we  denote  by  u^(x,t).  It  is  seen  immediately 
that  vij  satisfies 


3u^ 

Tt 


(>t  (  t  )  i  (  t  )  j  ^  (x ) 


(A-  2 4  ) 


The  perturbation  Wy  (t)  will  generally  oscillate  at  an 
acoustic  frequency  ,  so  Eq.  (A—  2  4)  predicts  that  u,  will 
oscillate  at  this  frequency  as  well.  This  mans  that A the  per¬ 
turbation  in  the  actual  field  E  will  oscillate  at  ,'0  f  \  .  ,  even 
though  tlie  cavity  mode  being  excited  would  oscillate  naturally 
(i.e.,  in  the  absence  of  the  perturbation)  at  the  frequency  . 

(since  it  is  degenerate  in  frequency  with  the  initially  present 
mode).  This  result  is  expected:  the  mode  N  is  in  a  state  of 
forced  oscillation,  and  it  oscillates  at  the  driven  frequency 

o.q  +  wA. 

Equation  (7e-24 )  is  in  the  form  of  the  st  irtitvi  equation  U  j 
the  flux  perturbation  in  our  treatment  in  iv.-et  ii  n  11.  it  is  t  In: 
seen  that  this  treatment  pertains  to  t  lie  flux  disturbance  assoc  Luted 
with  the  buildup  of  a  cavi  ty  mode  whose  uru-  :  ■  i  !":  >  ■  p.iency  i  s 

degenerate  with  that  of  the  initial  mode .  sisw  in  Section  Ii 

that  this  buildup  does  indeed  load  to  an  instability. 

On  the  other  hand,  we  now  consider  t  h<>  buildup  of  a  mode 
of  type  M,  satisfying  Kq .  (A-23).  Since  o.-ei  1  later,  at 

frequency  •.  A ,  Eg.  (A-23)  predicts  that  the“Lui  Idup  of  e.v  will 
be  resonantly  enhanced  when 


We  show  in  Sect  ion  IV  that  the  buildup  of  -•  node  c.f  this 
type  can  also  lead  to  an  instability,  at  least  in  the  low-fluenee 
limit. 
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APPENDIX  B 


PROPERTIES  OF 
FOR 


DEGENERATE  MODES  AND 
T11EIR  REIN FORC1 ; M E N T 


CONDITIONS 


We  wish  to  derive  several  results  stated  withovit  proof  in 
Section  II. 


First,  we  calculate  the'  conditions  under  which  two  resonator 
modes  may  be  degenerate  in  their  eigenf requeue ies .  Consider  for 
simplicity  light  modes  contained  in  a  box  of  length  L  and  width  W. 
The  eicjenf requrencies  (measured  in  cm-*)  satisfy 


(B-l  ) 


where  m  and  n  are  integers.  For  v  =  m/2L  >>  n/2VJ  we  may  write 
this  as  mo 


v 

mn 


v  +  ■= — 
mo  2  v 


mo 


(B-2) 


In  order  for  v  ,  .  to  equal  v  n  must  satisfy 

m- 1 , 0  ^  mn  J 


(B-3) 


(B-4) 


where  X  =  v  *  is  the  wavelength  of  the  light, 
mo 

Equation  B-4  suggest  r,  an  interr  n  on  t  h.  >  ord<  r  of  10  wi  1  1 
permit  the  frequency  "mn  of  a  t.ran  sverr.e  m<  to  b"  <i- cu  t  <■ 
with  that  of  a  longitudinal  mode.  For  oxamp  1  e ,  let  W  a  n:i, 

L  =  400  cm,  and  X  =  10"3  cm.  We  obtain 


n  *'  lG. 


(B-t) 
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From  this  value  wo  may  also  determine  the  number  n’  of 
transverse  oscillations  which  are  to  be  expected  in  an  acoustic 
wave  of  wave  vector  which  satisfies  the  phase-matching  con¬ 
dition 


k  ,  n 
m-1 , 0 


+ 


-► 


•  > 
k 

m ,  n 


(B-G) 


For  an  acoustic  mode  contained  in  the  transverse  dimension  of 
the  idealized  resonator  of  this  appendix  we  then  have 


v  -  'Ll  -  -fL7'  ’ll 

A  =  *A  =  W 

so  that 


(B-7) 


n ’  =  n/2  =  8.  (B-8) 

This  is,  of  course,  only  a  qualitative  result,  meant  to 
indicate  the  approximate  number  of  fringes  to  be  expected.  If 
the  approximate  conditions  of  the  model  (uniform  steady  state 
properties  for  flux  density  and  gain)  wore  rigorously  true,  tw 
would  have  to  be  a  integer  to  match  the  acoustic  boundary  condi¬ 
tions  at  the  walls.  In  the  actual  case  Eq.  (B-8)  is  only 
approx imute . 

Next  we  show  that  the  wave  scattered  by  the  phase  grating 
(standing  acoustic  wave)  produces  a  flux  perturbation  of  the 
same  shape  as  the  acoustic  wave. 

The  theory  of  scattering  of  3 ight  by  an  acoustic  wave 
(Debye-Scars )  predicts  scattered  or  diffracted  waves  at  angles 

of 


( B  -  9 ) 


with  respect  to  the  initial  wave.  These  waves.,  at  least  for  a 
confocal  unstable  resonator,  will  after  one  round  trip,  return 
at  the  prune  anjjl  e_s  . 


The  interference  between  those  wave;*,  and  the  unport  in  bed 
waves  will  be  of  the  form 


(B-10) 


if  uj  has-  the  same  import  urb-'d  f  r<  ; 1 i* -re "/  (  i  .  o.  ,  if  k]  k  )  as. 

flic  original  wave.  Th<*  s.  it  nation  is  l  1  1  u: .  1 1  . 1 1  •  ;  in  Figuii-  n- 1  . 
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cjurc  B- 1  .  Illustration  of  the  Peprodue  i  hi  1  i  ty  of  the  Scat  tered 
Wave.  The  unpot  t  urbed  v.mvc  focuses  (virtually)  at 
the  point  0  in  the  conCTun  f  oca  1  plant'  of  the  t  v;o 
mirrors.  The  sea*  torod  waver,  a  1  so  focus  in  this 
plane  hut  at  sliqht  ly  displa.cetl  points  A  and  A'  . 

The  return  waves  after  one  round  trip  w 1  1  ]  t  n.  },.* 
in  the  same  cl  i  root  i  on  as  the  or  i  Mina  1  scat  t  ered  wav*  r  . 
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Finally,  wo  justify  our  assertion  that  the  density  pertur¬ 
bation  is  independent  of  z,  the  light  propagation  direction.  If 
in  the  notation  of  Section  II  ui ,  and  hence  I]  ,  ha  5; ,  in  fact,  a 
different  k  ,|  than  ho,  then,  after  the  common  factor  exp  i  (kz  -  .ot 
is  factored  out,  4>1  J,  can  be  written 

*1£  ~  cos[kax  "  if]  (B“n) 

Now  the  gain  and  heating  equations  will  be  inconsistent  unless 
gj  i  and  p  j  p  have  this  same  additional  7.  behavior.  However,  sine- 
the  thickness  of  the  medium  is  of  the  order  of  half  the  separ.it'.  c: 
of  the  mirrors  (or  less)  in  a  typical  case,  the  /-variation  v.* ill 
be  insufficient  to  cancel  out  the  effects  of  the  mode-medium  into: 
action  and  we  are  justified  in  neglecting  it. 
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In  this  append  i  ,  we  present  a  semi-heuris ti  c  justifica¬ 
tion  for  Eq .  (60)  of  ! he  text. 

Let  us  assume  the  hen t inn  if;  oim  dimensional  (flow  direc¬ 
tion)  and  uniform  in  the  roe  ion  (-  b/2)  *  '  (b/2).  In  this 

case,  Fourier  and  Laplace  transforming  Eq .  (57)  loads  to 


(k ,  s) 


-  2(1  -  1)  F] (s)  k  sin  k  ~ 

(s  +  ivk )  [(s  +  ivk) 2  +  c  2k‘ 


(C-l  ) 


where  Fj(s)  represents  the  temporal  history  of  the  lieatjne  and 
wc  have  used  the  F-  transform  of  a  stir  function.  Then ,  by  con¬ 
tour  integration  we  obtain 


Pj  (x-0,s)  -  e 

C  S 


sb  _  ph  _ 

2v  2c.  'V  "”2 


sb  (  v/c 


cosh  -- 


(C-2  ) 


Expand  i  net  about  s  -  0  wo  find  tint  the  fired  two  terms  in  the 
brackets  are  the  same  as  the  expansion  of 


Pj  ( x  =  0  ,  s ] 


( Y  -  1  )  F  (  s. ) _ V1 

2  s  +  1 


(C-3) 


where 


1  -  -V-  - 

c  2 

s  c 


(C-'i) 


1^0 
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Ke  may  arrive  at  a  result  similar  to  L'q.  ( C  —  3  >  by  makina 
substitution 


s  +  ivk  -*•  s  +  v 2 


For  v/cg  ~  0.5,  we  have 


s  +  v_  x  s  +  1.25  p 
2  b 


In  the  text  wo  have  taken,  somewhat  arbitrary  s  +  v/b  as 
clompim?  term  on  intuitive  grounds  for  all  value's  of  v/cw.-, 
we  expect  the  convective  cl  am  pine  to  be  of  the;  order  of  ( i 

time) ~ * . 
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APPENDIX  D 


RESONATORS 


1.  EMPTY  RESONATOR  (STABLE.) 

There  are  many  ways  in  which  the  equations  satisfied  by 
the  eigenfunctions  of  an  empty,  stable  resonator  can  bo  formu¬ 
lated.  In  this  appendix  a  derivation  based  on  classical  poten¬ 
tial  theory  will  be  presented,  since  that  provides  the  most  lucid 
extension  to  the  analysis  to  be  presented  later  on. 

The  steady-state  wave  equation  for  a  transverse  electro¬ 
magnetic  wave  (Helmholtz  equation)  is  given  by  L2E  +  =  o 

where  k2  =  m2/c2 ,  is  the  wave  frequency ,  c  is  the  speed  of 
light  and  E  is  the  Laplace-transformed  transverse  component  of 
the  electric  field.  The  solution  to  this  equation  can  be  written 
in  terms  of  the  Huygens-Frcsnel-Kirchhof f  integral  by  means  of 
Green's  identities ( 1 ) 


E (x ,y , z) 


(D-l) 


where  G  is  the  appropriate  Green's  function,  A  is  the  area  of 
interest,  Es  and  DE  /Dn  are  the  electric  field  and  its  normal 
derivative  on  the  surface  respectively.  Thus,  once  the  field 
and  its  normal  derivatives  on  the  bounding  surfaces  are  given, 
the  field  everywhere  can  be  determined.  This  approach  lias  the 
disadvantage  chat  it  requires  the  spcci f i cat j on  of  both  the 
electric  field  and  its  normal  derivative  on  the  surface.  In 
general,  either  the  function  or  its  normal  derivative  will  be 
specified.  These  are  the  Dirichlet  and  Neumann  boundary  con¬ 
ditions,  respectively.  Note  that  mixed  conditions  are  allowed 
where  Dirichlet  is  specified  ever  part  of  the  bounding  surface 
and  Neumann  over  the  rest  of  it.  but  not.  both  at  the  same  time 
over  the  same  surface".  Thus,  it  is  necessary  to  reformulate 
the  solution  in  terms  of  so  called  dipole  and  singlet  charge 
distribution.  For  the  Dirichlet  problem  this  implies  that  (-) 


(1)  Born,  M.  and  Wolf,  E. ,  Principles  of  Optics,  Pergamon 
Press,  NY  (1959) . 

(2)  Mi  hill  in ,  S.G.,  Integral  Equations,  Peraamon  Press,,  NY 

(1957) . 
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A  /■ 


( x  ,  y  ,  z )  =  ^-  I  U  ( x  '  )  —  G(x,x')  dA 


(0-2) 


and  the  dipole  layer  ;i  satisfies 


E  —  + 


J"tn- 


™  dA 


(D-3) 


where  E  is  the  value  of  the  electric  field  on  the  boundary,  n 
is  the  normal  direction,  A  is  the  surface  area  and  G  is  the  ap¬ 
propriate  Green's  function  for  the  problem  at  hand.  For  three- 
dimensional  problems  it  is 


G  (x  ,  x  '  )  =  e 


ik  •  r 


(D-4  ) 


Consider  applying  this  procedure  to  the  stable  resonator 
shown  in  Figure  D-l  which  consists  of  the  two  mirrors  with  un¬ 
equal  radii  of  curvature  separated  by  a  distance  D.  In  the 
Fox- Li  calcualtions  for  the  stable  resonator  only  the  fields  on 
the  mirrors  are  determined,  that  is  only  n  is  solved  for.  Since 
fields  on  the  mirrors  vanish  Eq .  (D-3)  reduces  to 


‘l  +  j  ,J1  3n 


2  Dn . 


(D-5) 


u2  + 


f  Vl 


2  Sn. 


(D-G) 


where  U}  ,  vi 2  >  ,  M2  are  the  dipolo  layers  and  minor 

mirrors  1,  2,  respectively;  S}2»  S] \ ,  Sp2  nre  the  di 
tween  the  mirrors  including  curvature  turns  and  the 
of  different  points  on  the  same  mirror.  Writing  the 
terms  in  Hq .  (D-G)  ns  ,  1  ;i  1  and  the  first  and.  last  t 

(D-6)  as  Y2  H  2  these  oquations  can  be  rewritten  as 


a  rear 
st.  mo: 


of 


Vl 


'*  2  TTm 


(D-7) 
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Figure  D- 1.  Empty  Stable  Resonator 


avco  rvi  ni  it 


Substituting  Eq.  (D-8)  into  (D-7)  results  in  the  conventional 
integral  equations  for  the  Fox-Li  problem.  (3)  In  order  to  solve 
these  equations  the  paraxial  approximation  is  made,  that  is;,  the 
term  proportional  to  S~2  in  Eqs.  (0-7)  and  (D-8)  is  neglected. 

The  question  of  what  happens  when  the  resonator  is  loaded  will 
now  bo  addressed. 

2.  LOADED  RESONATOR  (STABLE) 

The  Max we 1 1  equations  appropriate  to  the  loaded  resonator 

are , 

V  •  e  E  =  -  V  •  P  ( D  -  9 ) 

o 


V  •  B  =  0  (D-10) 


V  X  E 


9B 

9t 


(0-11 ) 


V  x  H 


9E  , 
e  -vr  + 
o  6 1 


3  IE 

3t 


(D- 1 2 ) 


where  it  has  been  assumed  that  there  are  no  f roe-charge  currents 
flowing,  but  the  medium  responds  to  electromagnetic  waves  through 
the  polarization  term. 

Equations  (D-9)  and  (D-10)  can  be  combined  into  a  single 
wave  equation, 

,  2  2 

V2  E  -  — -  E  =  p  — — P  -  —  V  V  •  P  (D-13) 

2-2  O  -  .  2  t 

c  '  Dt  9t  o 


For  transverse  electromagnetic  waves  the  last  term  vanishes. 
The  polarization  term  is  assumed  to  be  given  by 


(3)  Fox,  A.G.,  Li,  T.  ,  Bell  Syst.  Tech.  J  40,  433  (191,1) 


I ^ 
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U5  =  c0  (x  ‘q;  +  c.c) 


(l)-]  •? ) 


where  x  is  the  complex  dielectric  susceptibility.  It;;  form  in 
terms  of  more  fundamental  quantities  is  given  in  A pg  ndix  ] 
Substituting  Eq.  { D  —  1 4 )  into  the  Laplace  transform  of  Kq .  ( 1 )  —  1  ) 

results  in 


V2  E  +  k2  E  =  0  ( D - 1 5 ) 


where 

k2  =  ~  (l  +  Xr  +  +  i  (D-1G) 

r  o 

where  p  denotes  a  phenomenological  loss  term  due  to  output 
coupling  from  the  resonator  ( ^  <  0  for  losing  medium)  .  The 
solution  to  Eq .  ( D— 1 6 )  can  be  written  in  the  same  form  as 
Eq.  (D-l).  However,  since  again  only  tire  field  or  the  normal 
derivative  would  normally  be  specified  it  is  necessary  to  re¬ 
cast  the  problem  in  terms  of  one  for  the  singlet  or  doublet 
charge  layers.  However  because  of  the;  amplification  by  the 
medium  it  is  now  necessary  to  calculate  the  fields  between  the 
mirrors  as  well.  The  only  difference  in  the-  doublet  equation 
(Eq.  ( 0 — 3 ] )  would  be  that  the  Green's  function  has  a  growing 
(decaying)  part  to  it  because  of  the  complex  values  of  k. 

Hence  special  attention  has  to  be  paid  to  the  convergence  prop¬ 
erties  of  the  various  functions  if  equations  of  the  form  of  Kqs. 
(D-7)  and  (D-8)  are  to  be  derived.  Most  of  the  difficulties  are 
swept  into  finding  the  appropriate  Green's  function. 

At  this  point  an  alternative  procedure  could  be  used. 

This  is  the  method  of  multiple  scaling.  This  technique  is 
used  to  essentially  eliminate  the  rapid  variation  in  t  he  /rN 
direction  (direction  normal  to  the  mirrors)  1  ir:n  tly  yqu.t  t  i  ons  . 

To  implement  this  the  fields  arc  assumed  to  have  e  varia¬ 
tions  in  the  'z'  direction.  Substitutin'!  this  form  into  Eq  .  (D-  1  5 ) 

and  averaging  over  the  fast  variation  results  in 

2 

V  2  E  +  2 ik  — -  E  +  k'2  E  ~  -  i  —  >  .  K  -  i,.  r  E  (D-17) 

T  z  dz  l  x  o  r 
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whore 


k 

7. 


2 


2 


2 


O' 


f 


(h- ] H) 


2 

is  the  transverse  Laplacinn  operator. 

The  solution  to  Eq .  (D-17)  can  in  principle  be  found  by 

the  following  procedure: 


E  (r) 


E  (  r  '  ) 


G  (V  ,  r  1  ) 


(0-19) 


(for  a  homogeneous  medium)  and  G  sotisfi 


G  +  2ik  G  -l  k  '  2G 

z  o  z 


(r-r  ’  ) 


(D-20 ) 


Note  that  E  and  G  have  to  be  found  in  t  ho  ent  i  re  r<  •:  ion  L>  -tva.  -n 
the  mirrors,  not  just  on  the  mirrors  themselves  as  had  been  t h*  • 
case  for  the  empty  resonator. 

Even  assuming  that  the  Green's  function  can  found 
Eq .  (D— 19)  is  still  an  integral  equal  ion  which  has  be 

solved,  either  numerically  or  in  a  perturb..- 1  :>  n  s  (i  cm 

appproximation )  .  Solving  the  coupled  sot.  of  E:s.  ('.-I'.1)  and 
(D-20)  is  more  akin  to  a  scattering  calculat  ion  the,  it  i 
using  tlie  Green's  theorem  of  the  previous  sect,  ie.n . 

3.  LOADED  RESONATOR  (UNSTABLE ) 


The  unstable  resonator  is  even  more  ces  s  lice-  i.  Eiuure 
D-2  illustrates  why  this  is  the  case.  The  Wuv •  if;.  (D-17)  '.till 
applies.  However,  it  has  to  be  aol  v<  d  in  a  num!  1  :  r-  :ics-- 
the  solutions  have  to  be  matched  across  t  h-  ■  bound  ;  :  s.  I  -1 )  - 

gion  I  contains  both  a  right-going  and  a  left  -goi  n-:  s.cse.  p.  ;n 
II  contains  only  a  r  i  g'ht.-goi  ng  wav--  with  cum.  ie-cion  III  e<  -n- 
tains  a  right  going  v.’av<>  but  with  no  sain  (.;  0).  in  t  1  • 
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where  f  and  G  arc  to  be  determined.  The  boundary  condition  E  =  0 
on  the  two  mirrors  lead  to 
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where  the  magnification  is  given  by  M  =  (D  +  d) /d  and  the  longi¬ 
tudinal  mode  number,  k  =  nr/D.  Substituting  Eq.  (D-21 )  into 
(D- 1 7 )  loads  to 
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whore  V»p2  is  the  transverse  I.aplacian  in  cylindrical  and  spher iea  1 
coordinates  in  the  two  equations  and  <j  -  -  (  .  2/k  ,.2cd }  k  >  i  is  tiio 
gain.  Equations  (0—24)  ,  and  (D-2r>)  were  solved  by  Moore  and 
McCarthy  (‘1 )  using  an  assumed  gain  profile.  The  self-consistent 
gain  profile  is  given  by  the  Rigrod  equation, 


V 


+  p 


(L>-2(. ) 


where?  V  is  the  flow  velocity  of  t  he  medium  ;  ,  is  the  deact  i  vat  ion 
time,  1  ,  s  are  the  optical  and  sat  ur.it  ion  flux  and  P  is  the  pump 
ing  rate .  Since  i  ~E'-  this  implies  that  the  oqu  i 1  i  enr.  lor  t  h*  • 
fields  [  (D—  24)  ,  ( U—  2  5  )  ]  are  coupled  thtouih  the  a  :  i  n  Eq.  (b-.M>). 

Since  the  two  fields  are  do.-.cr  ibed  in  d  i  1  1  •  u  t  :i  t  c  -u  i :  n  ■  t  e  r.y.;- 
toms  ,  this,  mala  s  the  se  1  f -cons  i : .  t  <  n  t  pm!  ■ !  *  m  put  lcnlni  1  y 
plcx.  A  complete  numerical  analys.  i s.  of  t  lit -se  Kgs.  (!)-.!<» )  ,  (P-.>:t) 

and  (D-26)  is  described  in  Sect,  ion  V. 
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In  the  integral  formulation  of  the'  mode-medium  inntubi  1  i  t 
the  solution  to  the  equilibrium  problem  is  written  as  uQ  :  /K(  u 
where  Kq  is  the  appropriate  Huygens-Fresnel  operator. 
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APPENDIX  E 


DERIVATION  OF  THE  POLARIZATION'  TERM 


In  the  simple  approach  taken  here  the  electrons  are  treat 
somiclasy  .ical  ly .  Since  the  medium  treated  is  a  1  ow-dens  i  ty  gas 
the  electric  field  seen  by  the  individual  molecules  is  the  same 
as  the  external  field.  (In  a  denser  medium  like  a  solid  the  ei 
tive  field  would  be  Ecff  “  E  +  IT/ 3  where  IP  i  s  the  polarization 
shielding  .  )  ( 1 ) 

Thus  individual  electrons  obey  the  following  equal  ion , 

X  4  ,X  a  2  X  -  ;;  E  (E- 


where  y,  3  arc  the  damp  i  no  (rrowth)  rat  e,  the  orbit  a  1  bin  dir  1 
frequency  and  E  is  the  s  ti  mu  1  a  t  i  ng  eb-etiic  field.  Solving  Kg, 

(E- 1 )  leads  to 

-  i  .<  t 
E.  e  ° 

x=  2--~  --2— n:- 

L  ‘  o  B  J  1  ‘  o  -* 

where  c-c  is  the  complex  conjugate.  The  dipole  moment  of  the 
molecule  is  given  by 

p  =  -ex  (E- 

Finally  the  macroscopic  polar’  ation  for  a  density  of  N  mol  ecu 
per  cubic  centimeter  is  give. 
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The  effects  of  a  spatial  ly  vaiyinq  wave  are  i  in' !  1;  ! 
fashion.  Th>>  Fourier  transformed  polar  i  zat  i  i<n  he- 
nary  parts  given  by, 


ic  !  1:  !•  ‘  in  a  s  i m 


1  1  an:  1  i . 
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where  R  =  ( u>02  -  !b  2 )  2  +  -y  20.o2 .  Thus  the  real  part  modifies  the 
index  of  refraction  while  the  imaginary  part  produces  st  imul  atod 
emission  of  absorption  depending  upon  the  population  inversion. 

In  Section  IV  IP  is  written  as  follows: 

IP  =  c Q x r E  +  ieQ  xA  E  (Id- 7 ) 


Comparing  Eq.  ( E—  7 )  with  Eqs.  ( E—  5 )  and  (E-6)  leads  to 
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(Id- 9) 


In  addition  IP  was  defined  in  terms  of  a .  Equations-  (E-S)  and 
(E—  9)  give  a  simple  definition  of  the  frequency  dependence  of  t . 


In  the  semic  1  assica  1  treatment  growth  (stimulated  emission) 
is  obtained  by  lotting  r  be  negative. (2)  In  addition,  the  effect: 
of  saturation  enter  y  by  making  it  amplitude  dependent.  Thus  a 
relationship  between  the  sem  iclass  ica  1  treatment,  .uvi  the  quantum 
treatment  is  obtained  by  identifying  , ;j  with  the  gain  of  the 
medium.  The  gam  is  given  by 
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Hence 
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Finally  this  allows  the  p<d  u nwtion  to  be  written  in  trims  oi 
the  gains  as. 
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(K-12) 
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This  result  is  derived  on  the  basis  of  assumin'-!  that  M  and  q 
remained  constant  over  the  time  period  of  interest.  The  generali¬ 
zation  allows  both  of  these  to  change  and  is  treated  in  Section 
IV.  Note  that  even  in  equilibrium  g  is  allowed  to  vary  because 
of  saturation  dependence  on  the  fields.  That  is  y  depends  on  t. ho 
fields,  but  is  assumed  to  have  a  weak  dependence  on  the  fields 
at  Lest.  For  single  frequency  electric  fields  the  polarization 
in  real  space  and  time  is  given  by 

IP(rf  t)  =  IPk  +  IP  *  (H-l  3) 

where  IP^ is  given  by  Eq.  (E-12).  For  more  complicated  electric 
fields  the  Fourior-Laplace  inversion  of  Eq .  (E-12)  has  to  be  per¬ 

formed.  This  leads  to  a  convolution  integral  of  the  form, 

IP  (r  ,  t)  =  [dr  ’dt1  x(r’,t')  E(r-r’,t-t')  (E-14) 

where  X  is  given  by  Eqs.  ( E—  8 )  and  (E-9). 
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